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PREFACE 


The object of this book is to provide a compact ex¬ 
position of the fundamental results in the theory of 
tensors and also to illustrate the power of the tensor 
technique by applications to differential geometry, 
elasticity, and relativity. In the first five chapters the 
mathematical concepts are developed without undue 
stress on rigour. The remaining three chapters are in¬ 
dependent of one another except that sections 88 and 89 
of chapter VI, which treats Euclidean three-dimensional 
differential geometry, are necessary for a proper under¬ 
standing of chapter VII which contains the theory of 
cartesian tensors and elasticity. Finally, chapter VIII is 
devoted both to the special and general theories of 
relativity. In the limited space available it is impossible 
to do justice to the physical principles underlying both 
these theories. But in order to help the reader unac¬ 
quainted with relativity some explanatory matter has 
been incorporated into the text. 

The presentation owes much to the authors listed in the 
bibliography, especially to McConnell, Synge and Schild. 
In particular, I wish to express my thanks to Dr. D. E. 
Rutherford for numerous suggestions and helpful 
criticisms during the manuscript and proof stages. Lastly 
I wish to thank my wife for her help with the proof¬ 
reading. 

B. S. 

Trinity College, Dublin; July, 1952 

PREFACE TO THE THIRD EDITION 

In this edition various errors have been corrected. 
Further, I wish to thank Mr. L. Lovitch for a neat proof 
(inserted on page 56) that the curvature tensor is zero in 
a flat space. 


B. S. 
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A Concise Course 



CHAPTER I 


TENSOR ALGEBRA 

§ 1. Introduction 

The concept of a tensor has its origin in the develop* 
ments of differential geometry by Gauss, Riemann and 
Christoffel. The emergence of Tensor Calculus, otherwise 
known as the Absolute Differential Calculus, as a syste¬ 
matic branch of Mathematics is due to Ricci and his pupil 
Levi-Civita. In collaboration they published the first 
memoir on this subject: — * Methodes de calcul differential 
absolu et leurs applications', Mathematische Annalen, vol. 
54, (1901). 

The investigation of relations which remain valid when 
we change from one coordinate system to any other, is 
the chief aim of Tensor Calculus. The laws of Physics 
can not depend on the frame of reference which the phy¬ 
sicist chooses for the purpose of description. Accordingly it 
is aesthetically desirable and often convenient to utilise 
the Tensor Calculus as the mathematical background in 
which such laws can be formulated. In particular, 
Einstein found it an excellent tool for the presentation of 
his General Relativity theory. As a result the Tensor 
Calculus came into great prominence and is now invalu¬ 
able in its applications to most branches of Theoretical 
Physics; it is also indispensable in the differential geome¬ 
try of hyperspace. 

It is assumed that the reader has an elementary know¬ 
ledge of determinants and matrices. As he may not be 
acquainted with the Calculus of Variations, the minimum 
problem in the theory of geodesics is treated from first 
principles. 
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§ 2. N-Dimensional space 

Consider an ordered set of N real variables x\ x\ ..., 
x { ,..x N ; these variables will be called the coordinates 
of a point. (The suffixes 1 , 2, ... i, ... N, which we shall 
call superscripts, merely serve as labels and do not possess 
any significance as power indices. Later we shall introduce 
quantities of the type a t and again the i, which we shall 
call a subscript, will act only as a label.) Then all the 
points corresponding to all values of the coordinates are 
said to form an N-dimensional space, denoted by Vs- 
Several or all of the coordinates may be restricted in 
range to ensure a one-one correspondence between points 
of the Vs and sets of coordinates. 

A curve in the Vs is defined as the assemblage of 
points which satisfy the N equations 

x* = x*{u), (i = 1, 2 , -. - N) 
where tt is a parameter and x*(u) are N functions of tt, 
which obey certain continuity conditions. In general, it 
will be sufficient that derivatives exist up to any order 
required. 

A subspace Vm of Vs is defined for M < N as the 
collection of points which satisfy the N equations 
x 1 = x^u 1 , tt 8 , . . . u u ), (i = 1, 2,... N) 
where there are M parameters tt 1 , « 3 ,... u u . The 
u 8 ,... u M ) are N functions of the tt 1 , tt 8 ,... u u 
satisfying certain conditions of continuity. In addition 
the M x N matrix formed from the partial derivatives 
dx^du 1 is assumed to be of rank M*. When M = N — 1, 
the subspace is called a hypersurface. 

§ 8. Transformation of coordinates 

Let us consider a space Vs with the coordinate system 
x\ x 2 ,... a^. The N equations 
(8.1) x 1 = 9 ? < (® 1 , a 8 ,... (t = 1, 2, . .. N) 

* T. Levi-Civita, The Absolute Differential Calculus, pp. 0—12. 
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where the q>* are single-valued continuous differentiable 
functions of the coordinates, define a new coordinate 
system x 1 , x 8 ,... x^. Equations (8.1) are said to define 
a transformation of coordinates. It is essential that 
the N functions q> { be independent. A necessary and suf¬ 
ficient condition is that the Jacobian determinant formed 
from the partial derivatives dx i fdx i does not vanish*. 
Under this condition we can solve equations (8.1) for the 
x* as functions of the x* and obtain 

x* = ^(x 1 , S 2 ,... 5^) (i = 1, 2, ... N ). 


§ 4. Indicial and summation conventions 
We will now introduce the following two conventions: 

(1) . Latin indices, used either as subscripts or superscripts, 
will take all values from 1 to N unless the contrary is 
specified. Thus equations (8.1) are briefly written 
x* = ^(x 1 , x 2 , ... aF), the convention informing us that 
there are N equations. 

(2) . If a Latin index is repeated in a term, then it is un¬ 
derstood that a summation with respect to that index over 
the range 1, 2,... N is implied. Thus instead of the ex- 

N 

pression Za { x*, we merely write a i x i . 

<=i 

Now differentiation of (8.1) yields 

N P) m { N 

dx* = Z ^- r daf = Z t-daf, (i *> 1 , 2,. .. N) 

r=l dx T r=l ax T 

which simplify, when the above conventions are used, to 

(4.1) dx i = ^- r dx r . 

The repeated index r is called a dummy index, as it 
can be replaced by any other Latin index, except i in 

* R. P. Gillespie* Partial Differentiation* pp. 48—46. 
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this particular case. That is, equations (4.1) can equally 

well be written dx* = • 5 — dx m or for that matter 

dx 1 — •=—; dx 1 . In order to avoid confusion, the same 
dx* 

index must not be used more than twice in any single 

N 

term. For example (Z a { x*) a will not be written a^a^* 

<~i 

but rather a i a j x i x*. It will always be dear from the con¬ 
text whether x 2 means x with superscript 2 or a squared. 
Usually powers will be indicated by the use of brackets; 
thus (x 1 *) 2 means the square of aF. The reason for using 
superscripts and subscripts will be indicated in due 
course. 

Let us introduce the Kronecker delta defined by 


(4.2) 


I sf=i aj=k 
l if = o if j # *. 

An obvious property of the Kronecker delta is that 
djA 1 = A k , since in the left-hand side of this equation 
the only surviving term is that for which j — ft. Also 
dx k [dx i = 6f, because the coordinates x* are independent. 


Ex. Show that <5^5^ 


4i 6 \ = N ‘> 


'b&'b£ { 

IS 1 


§ 5. Contravarlant vectors 
A set of N functions A* of the N coordinates x i are said 
to be the components of a contravarlant vector if they 
transform according to the equation. 

< 6 -» 

on change of the coordinates x* to x*. This means that 
any N functions can be chosen as the components of a 
contravariant vector in the coordinate system X 1 , and the 
equations (5.1) define the N components in the new 
coordinate system x*. On multiplying equations (5.1) by 
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dap fix 1 and summing over the index i from 1 to AT, we obtain 
dx* t. dx k dx 1 dap . k Ale 

W A= U‘W a, = 1 >S a =i ' A=A - 

Hence the solution of equations (5.1) is 

(5-2) A '=% A ‘- 


When we examine equations (4.1) we see that the 
differentials dx 1 form the components of a contravariant 
vector, whose components in any other system are the 
differentials dx 1 of that system. It follows immediately 
that dx^du is also a contravariant vector, called the 
tangent vector to the curve x { == x*{u). 

Consider now a further change of coordinates 
x'* = g^x 1 , x 8 ,... x"). Then the new components A H 
must be given by 


A'* = — Ai = 
dxP 


d J?l ^ A* = — A* 

dx-' dx* daP ' 


This equation is of the same form as (5.1), which shows 
that the transformations of contravariant vectors form 
a group*. 

With the exception of the coordinates x l themselves 
a single superscript will always denote a contravariant 
vector unless the contrary is explicitly stated. The coor¬ 
dinates x i will only behave like the components of a 
contravariant vector with respect to linear transfor¬ 
mations of the type £* = where the a\ are a set 
of N 2 constants, which do not necessarily form the com¬ 
ponents of the entity introduced in section 8 and there 
called a tensor. For in this case d&jdx* = a] and the trans- 

formation can be rewritten x { = xK With respect to 

ox* 

general transformations of coordinates, the x 4 do not form 


* W. Ledermann, Introduction to the Theory of Finite Groups, 

pp. 2—8. 
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the components of a contravariant vector. Essentially 
this means that if we select A { = x*, then the new com¬ 
ponents A* with respect to the coordinate system x* do 
not satisfy the equations A* = x { . 


Ex, If a vector has components 


d*x d t y 
dP* di * 


in rectangular cartesian 


coordinates, show that they are ——r — + — — — in 

\dtl <U* ^ r di di 

polar coordinates. 


§ 6 . Covariant vectors 

A set of N functions A i of the N coordinates x i are said 
to be the components of a covariant vector if they 
transform according to the equation 


( 6 . 1 ) 



on change of the coordinates x* to 5 *. Any N functions 
can be chosen as the components of a covariant vector 
in the coordinate system x\ and the equations ( 6 . 1 ) 
define the N components in the new coordinate system 
On multiplying ( 6 . 1 ) by dx { fdx k and summing over the 
index i from 1 to N, we obtain 


( 6 . 2 ) 


dx l j _ dx* dx? A doc? A 

dx k Ai ~dx* w*'* ~~W kA *~ Ak ' 


0 / 0 / 0 ®* . 

Smee ^ ^ 7 , it follows immediately from ( 6 . 1 ) 


that the quantities 0 // 0® 1 are the components of a co- 
variant vector, whose components in any other system 
are the corresponding partial derivatives df/dxK Such a 
covariant vector is called the gradient of /. 

A single subscript will always denote a covariant vector 
unless tiie contrary is explicitly stated. In conformity 
with this convention we shall regard the index i in the 
covariant vector djjdx 1 as a subscript. 
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57 


We now show that there is no distinction between con- 
travariant and covariant vectors when we restrict our¬ 
selves to transformations of the type 


(6.3) x ( = a i m x m + b*, 

where b* are N constants which do not necessarily form 
the components of a contravariant vector and a { m are 
constants (not necessarily forming a tensor) such that 

« = 

We multiply equations (6.8) by a*, and sum over the 
index i from 1 to N and obtain 


Thus 


x T — a^x* — a^b*. 
dx l _ dx* _ 

d^~W { ~ ah 


which shows that the equations (5.1) and (6.1) define 
the same type of entity. 

Ex. Prove that the transformations of covariant vectors form a 
group. 


§ 7. Invariants 

Any function I of the N coordinates x* is called an 
invariant or a scalar with respect to coordinate trans¬ 
formations if 1=1, where 1 is the value of I in the new 
coordinate system x { . 

From the components A* and B { of a contravariant and 
covariant vector respectively, we can form the sum A i B i . 
When we change to new coordinates x*, this sum trans¬ 
forms to A i B i . Now 

ZE, = §§ * gj£ B, = 9lA>B t = A*B t . 

That is, 

AiBi = A*B { . 

Thus A*B { is an invariant. 
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Another invariant is the quantity 


^2 • • • 4* <5$ = AT. 


§ 8 


§ 8. Second order tensors 

Form the N 2 quantities A if = BW, where B* and C i 
are the components of two contravariant vectors. It 
follows from (5.1) that the A i} transform according to 


( 8 . 1 ) 


A*> = 


A* 

daP da j* 


More generally, if we have N 2 functions A ij whose 
transformation law is that of (8.1), then we call A** the 
components of a contravariant tensor of the second 
order. This tensor is not necessarily the product of the 
components of two contravariant vectors. Any set of N 2 
functions can be chosen as the components of a contra¬ 
variant tensor of the second order, and then (8.1) defines 
the components in any other coordinate system x l . 

Similarly, if we have N 2 functions A ti whose transfor¬ 
mation law is 


( 8 . 2 ) 


j _ dx k da? . 

A{i ~ W< d& A * 1 ’ 


we call A {j the components of a covariant tensor of the 
second order. 

Further, if we have N 2 functions Aj whose transfor¬ 
mation law is 


( 8 . 8 ) 


4 - 


dx { dx l k 
dx*d& Al ’ 


we call Aj the components of a mixed tensor of the 
second order. 

Note that the indices are placed on the tensors as super¬ 
scripts when they denote contravariance and as subscripts 
when they denote covariance. In particular, the mixed 
tensor Aj transforms like a contravariant vector with 
respect to the index i and like a covariant vector with 
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respect to the index j. Consequently the i is placed as a 
superscript whilst the j is placed as a subscript. 

Now choose A\ to be the Kronecker delta df. From 
(8.8) we have 

Zi — _ dx { dx k _ dx* _ 

* dx k dx* 1 dx k dx* dx s * ’ 

That is, the Kronecker delta is a mixed tensor of the 
second order whose components in any other system 
again form the Kronecker delta. This justifies the 
placing of one index as a subscript and the other as a 
superscript. Yet if we select the N a quantities d if — <3J 
as the components of a covariant tensor in a coordinate 
system x*, the components in the new system x* are 

given by d if = -^r. . Therefore the transformed com- 
dx* dx’ 

ponents 6 {i do not form the Kronecker delta. 

Ex. Prove that A ij B , C s is an invariant, if B f and C* are contra- 
variant vectors and A ti a covariant tensor. 


§ 9. Higher order tensors 
A set of N *+* functions of the N coordinates 

x* are said to be the components of a mixed tensor 
of the (s -f- p)-th order, contravariant of the s-th order 
and covariant of the p-th order, if they transform ac¬ 
cording to the equation 


(9.1) 



dx u i dx u * dx 9 1 daftr 
dx t i " ’ " dx** dx T i * * * dx T v 


on change of the coordinates x { to x*. This formula, 
although rather formidable in appearance, is merely a 
combination of (5.1) with respect to contravariant in¬ 
dices and of (6.1) with respect to covariant indices. 

The order of indices in a tensor is important. The 
tensor A** is not necessarily the same as the tensor A t{ . 
(In the language of matrices A ){ is the transpose of A* 1 .) 
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If however two contravariant indices or two covariant 
indices can be interchanged without altering the tensor, 
it is said to be symmetric with respect to these indices. 
We shall now prove that if a tensor is symmetric with 
respect to two indices in any coordinate system, it 
remains symmetric with respect to these two indices in 
any other coordinate system. There is no loss in generality 
in proving this for the contravariant tensor A ij = A li . 
On applying (8.1) we have 


<»•*> 


dx s dx* 
dx l dx k 


A » = A 1 *, 


as required. We cannot usually define symmetry with 
respect to two indices, of which one denotes contra- 
variance and the other covariance, because this symmetry 
may not be preserved after a coordinate transformation. 
The Kronecker delta, however, is a mixed tensor which 
possesses symmetry with respect to its two indices. 

When all the indices of either a contravariant or a 
covariant tensor can be interchanged without altering 
the tensor, it is said to be symmetric. A symmetric tensor 
of the second order has at most $N(N + 1) different com¬ 
ponents. 

A tensor each component of which alters in sign but not 
in magnitude when two contravariant indices or two co¬ 
variant indices are interchanged is said to be skew- 
symmetric with respect to these indices. It can be 
shown by equations similar to (9.2) that the property 
of skew-symmetry is also independent of the choice of 
the coordinate system. Skew-symmetry, like symmetry, 
cannot be defined with respect to two indices, of which 
one denotes contravariance and the other covariance. 

If all the indices of a contravariant or a covariant 
tensor can be interchanged so that the tensor changes its 
sign at each interchange of a pair of indices, the tensor 
is said to be skew-symmetric. A skew-symmetric tensor 
A lt of the second order has at most \N{N — 1) different 
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arithmetical components, as all the quantities A u (no sum¬ 
mation) are zero. The several components of a skew- 
symmetric tensor of the N-th order are either zero or differ 
merely in sign. So there is essentially only one non-zero 
component of such a tensor. 

The most important deduction from (9.1) is this: if all 
the components of a tensor in one coordinate system are 
zero at a point, they are all zero at this point in every 
coordinate system. Further, if the components are iden¬ 
tically zero in one coordinate system, they are also iden¬ 
tically zero in every coordinate system. It is this property 
which constitutes the importance of tensors in physical 
applications. 

When a tensor is defined at all points of a curve or 
throughout the space V N itself, we say that it constitutes 
a tensor-field. 

Ex. 1. If A u is s skew-symmetric tensor, prove that 
(<5<<5? + 6\6*)A tb = 0. 

Ex. 2. Prove that the transformations of tensors form a group. 

§10. Addition, subtraction and multiplication of 
tensors 

It is clear that we cannot expect to give any tensorial 
meaning to the expression A u -j- B* t because it cannot 
satisfy the transformation law (9.1). It does, however, 
follow from this equation that any linear combination of 
tensors of the same type whose coefficients are invariants 
is a tensor of the same type. For example, from the two 
tensors Aj k and Bj k , we can form the tensor XA\ k + fiBj k 
which will satisfy (9.1) provided that A and ft are in¬ 
variants. In particular, Aj k -f Bj k and A$ k — Bj k are 
called the sum and difference respectively of the two 
tensors. As another example, we can write 

■Au = i(^« + A j{ ) + ^(Aij — A it ). 

Now A ti + A H is symmetric and A it — A it is skew-sym- 
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metric. Thus any covariant tensor of the second order 
is the sum of a symmetric and a skew-symmetric tensor. 
This is, of course, also true of a contravariant tensor of 
the second order. 

Let us select two tensors, one of contravariant order 
a and covariant order p, the other of contravariant 
order t and covariant order q. It then follows from (9.1) 
that the products of the components form a mixed tensor 
of contravariant order s + t and covariant order p + q. 
This tensor is called the outer product of the two 
tensors. For instance A*,^ nt = B% C^ n< is the outer pro¬ 
duct of the two tensors B*J and C l mnt and is a tensor of 
the type indicated by its indices. 

The division, in the usual sense, of one tensor by another 
is not defined. 

§11. Contraction 

Let us start with any mixed tensor, say A^ mn> and form 
the sum A^ mi . From (9.1) we have 

jgt __ dx* dx 1 daf dap* daP 0 
pvr das* dx* dx p dx Q dx T tmn ’ 

Therefore 

Tjr _ 9®* dx T dot daif* daf 1 {f 
Apv ~fa i dxid&d&dv A ' mn 
= d&dt*_ fa* u 
dx* dx p dx q i lmn 
_ d& daf dx m {n 
dx* dxP dx q lmn ' 

Thus we observe that A\ ! mj is a mixed tensor, contra¬ 
variant of the first order and covariant of the second order. 
This process, which is called contraction, enables us to 
obtain a tensor of order r — 2 from a mixed tensor of 
order r. In the above example we could contract a stage 
further and arrive at the covariant vector Afy. When 
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contracting, any superscript may be used to sum with 
any subscript. Therefore we can form the following dif¬ 
ferent tensors by contraction:- A^ mj , A^, Aj^, A^ mV 

Ail Ail Ai! Ail Ail Ait Ail Ail Tf 

■“»n» “(mu’ ■“!«» ■“«<> ^iml* A tmi* ■“iln “ uu ^lin' “ 

the tensor A" mn possesses any symmetric properties, 
there will be fewer tensors formed from it by con¬ 
traction. As another example, the invariant Aj is formed 
by contraction from the mixed tensor Aj. Tins justifies 
us in calling an invariant a tensor of zero order. 

We can also combine multiplication and contraction 
to produce new tensors. From the tensors A% and B l mnt , 
we may obtain such tensors as A*J A% B\ n} , Aff B^ j{ 
and many others. This process is called inner multipli¬ 
cation of two tensors and the resulting tensor is called 
an inner product of the two tensors. 

Note carefully that we never contract two indices of 
the same type as the resulting sum is not necessarily a 
tensor. Also it should now be clear that with our index 
notation, the summation convention generally applies 
to two indices one of which is a superscript and the other 
a subscript. 

§12. Quotient law 

Sometimes it is necessary to ascertain whether a set 
of functions form the components of a tensor. The direct 
method requires us to find out if they satisfy a tensor 
transformation equation of the type (9.1). In practice 
this is troublesome and a simpler test is provided by the 
quotient law. The quotient law states that N* functions 
of x 1 form the components of a tensor of order p, (whose 
contravariant and covariant character can readily be 
determined), provided that an inner product of these 
functions with an arbitrary tensor is itself a tensor. It 
will suffice to set out the proof for the following parti¬ 
cular case. The set of N 3 functions A i,k form the com¬ 
ponents of a tensor of the type indicated by its indices if 
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A” k B> = C*\ 

provided that is an arbitrary tensor and C pk a tensor. 
The transformed quantities, referred to a system of coor¬ 
dinates x*, satisfy the equations 

= C**, 

which become on substitution from (9.1) 


x m diP dxm _ dz> dx* 

dx l dx* dx* mn dcfi da? 

= diP dx k jiff.jp 
da?da? A 


With a change of dummy indices, we have 


dx» 

dx 1 



dx* dx* dx T \ 


K. 


= 0 . 


On multiplying this equation by da?jdxP and summing 
over p from 1 to N, (in future we shall merely write ‘on 
inner multiplication by dx’jdx*') we obtain 




Kn = 0. 


Since B* mn is an arbitrary tensor, we can arrange that 
only one of its components differs from zero. Now each 
component of B‘ mn may be selected in turn as that one 
which does not vanish. This shows that the expression 
in brackets is identically zero. That is 


( 12 . 2 ) 


A** 


dx™ dx n _ .. | | r dx k 
di*d&~ A da? 


Inner multiplication of this equation by ^ yields 
the result 


3x* dx* d& 
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Thus A mnr is a tensor of the third order and contra- 
variant in all its indices. In the above proof it is important 
that the tensor B l mn shall be arbitrary and must not 
possess any symmetric or skew-symmetric properties. Let 
us examine what happens if B l mn is symmetrical in m 
and n. It is no longer possible to deduce (12.2) from (12.1). 
The correct deduction now is that 


_ . 0®"* dot* r ... dx n daF* . 35* . 

Jiik -L. liik _ u _ Amnr IZ_ _]_ Anmr _ 

dx { dx* ' Z»r ^ a.-rr * 


dx* d& 


daf 


dx T 


On changing several dummy indices, this equation be¬ 
comes 




Inner multiplication by ^ shows that A mnr + A nmr 

is a contravariant tensor of the third order. If in 
addition we know that A mnr is symmetric in m and n, 
it follows immediately that A mnr is a symmetric tensor 
with respect to the indices m and n. We observe, there¬ 
fore, that the quotient law must be applied with care. 

Ex. 1. If A* and B* are arbitrary contravariant vectors and 
C ti A i B l is an invariant, show that C w is a covariant tensor 
of the second order. 

Ex. 2. If A* is an arbitrary contravariant vector and C it A t A l is 
an invariant, show that Cu + Ch is a covariant tensor of 
the second order. 


§ 18. Conjugate symmetric tensors of the second 
order 

Consider a symmetric covariant tensor of the second 
order A ti whose determinant | A t) | ^ 0. Let B (i denote 
the expression formed by dividing the cofactor of A it in 
the determinant | A ti | by | A u | itself. We shall prove 
that the B il so obtained are the components of a contra¬ 
variant tensor of the second order. In anticipation of 
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this, B fi is labelled as if it were a contravariant tensor. 
We have from the theory of determinants*. 

(13.1) A if B ik = dt. 

We cannot establish the tensor character of B ik by 
applying the quotient law directly to this equation, 
because A if is not arbitrary. Let us choose an arbitrary 
contravariant vector C f . Then D t = A (i C f is an arbitrary 
covariant vector, because theseNequations can be unique¬ 
ly solved for the C* in terms of the D { since | A i} | ^ 0. 
Consequently, 

D t B ik = A if a B {t = 8$ C> = CK 

Now if we apply the quotient law to the equation 
D t B ik = C k , we see that B ik is a contravariant tensor of 
the second order. Also it is clear from the definition that 
B u , like A ijy is symmetric. 

We will next attempt by the same process to obtain 
another tensor from B ik . Let E {f denote the cofactor 
of B li in the determinant | B ti | divided by | B u | itself. 
From the theory of determinants | A u |. | B a | = 1 and 
consequently | B t) \ =£ 0, which means that E (i always 
exists. Further we have 

E f ,B ik = <5?. 

Inner multiplication by A kl yields on application of 

(18.1) that 

Eif = Afi = Ai t . 

Thus this process only leads back to the original covariant 
tensor of the second order. We say that A if and B li are 
conjugate tensors if they satisfy equations (18.1). It is 
important to note that a tensor of the second order has 
a conjugate only if its determinant is not zero. 

Ex. If A (i = 0 for i ^j, show that the conjugate tensor B‘> = 0 
for t and B li = 1}A U (no summation). 


* A. C. Aitken, Determinants and Matrices, pp. 51—52. 
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§ 14. Fundamental tensor 
At this stage we introduce the concept of distance into 
our space V^. If the distance ds between the neigh¬ 
bouring points with coordinates x* and x* + dx { is given 
by the quadratic differential form 

(14.1) ds* = g ii dx { dx i 

where the g u are functions of x l , subject only to the 
restriction g = \ g u | ^ 0, the space is said to be a 
Rlemannlan space. In addition we postulate that the 
distance between two neighbouring points is independent 
of the coordinate system. That is, ds is an invariant. 
From the quotient law, it follows that g u + in is a co¬ 
variant tensor of the second order. We can write 

8n — h(in + in) + £(g« — in)- 

The contribution of $(g (i —• g ji )dx i dx i to ds* is zero, hence 
there is no loss of generality in assuming that g fl is sym¬ 
metric. Thus g if is a covariant symmetrical tensor of the 
second order called the fundamental tensor of the Rie- 
mannian space. The quadratic form g ij dx i dx i is called the 
metric. It is also the square of the line-element ds. 

The line-element ds of a three dimensional Euclidean 
space, referred to a system of rectangular cartesian axes, 
is 

ds* = (d* 1 ) 2 -f (d* 2 ) 2 + (d* 3 ) 2 . 

All the components of the fundamental tensor are zero 
except g u = gjg = g*j 3 = 1. It is evident that the metric 
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of a Euclidean space is positive-definite*. That is, da 2 
is zero when d® 1 = d® 3 = d® 3 = 0 but can only take 
positive values for all other real values of dx l , d® 3 and d® 3 . 

The Special Theory of Relativity discusses the four 
dimensional space with the line-element da given by 

(14.2) da 3 = — (d® 1 ) 8 — (d® 8 ) 8 — (d® 8 ) 8 + cW) 8 

This metric is not positive-definite, as it is positive for 
all curves along which x\ ® a and ®* are all constants, but 
negative for all curves along which ®* is constant. Thus 
along these latter curves the distance between neigh¬ 
bouring points cannot be real. In order that the distance 
da between two neighbouring points be real, equation 
(14.1) will be amended to 

(14.8) ds a = eg {i dx l dx s , 

where the factor e called the indicator takes the value 
-f 1 or — 1 so that da 2 is always positive. 

Ex. Show that the metric of a Euclidean space, referred to spherical 
polar coordinates as* = r, x* — 0 and *» = y> is given by 
ds* s dr* + r 3 dd 2 + r* sin® 0dy>®. 

§ 15. Length of a curve 
Consider the curve x* = x*(t) with the parameter t. 
From (14.8) the length of the curve between the points 
corresponding to t = ty and t = t% is given by 

(“) —\ h l «»*■**• 

If ^ ^ = o along a curve, then the two points 

corresponding to ty and are at zero distance from one 
another, despite the fact that they are not coincident. 


* A. J. Me. Connell, Absolute Differential Calculus, p. 10. 
A. C. Aitken, Determinants and Matrices, p. 187. 
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Such a curve is called minimal or null. The curve 
given by 


(15.2) 


x 1 = cjr cos d cos ip dt, at 2 = cj r cos 0 sin ip dt, 
a? = cjrswOdt, ar 1 = Jrdi, 


where r, 6 and ip are functions of t is a real null curve 
in the F 4 whose metric is (14.2). It is clear that no real 
null curves lie in a Riemannian space whose metric is 
positive-definite. 

A curve will consist of portions along which the in¬ 
dicator e is + 1, portions along which the indicator is 
— 1, and null portions. The length of the curve is then 
the sum of the lengths of these portions, the null part 
contributing zero to the value of the length. 

Except in the case of null curves, the parameter t may 
be chosen as the arc-distance a from some fixed point 
of the curve. From (14.8) 


(15.8) 


dx i dx* 

g {t ds da ~ 6 


along any portion of a curve which is not null. 


§16. Magnitude of a vector 

The magnitude A of the contravariant vector A* is 
defined by 

(16.1) (A)* = e {A) g fi A*A*, 

where e^) is the indicator + 1 or — 1 which makes A 
real. The magnitude A is an invariant. In a Euclidean 
space, referred to rectangular cartesian coordinates, 

(16.1) reduces to the familiar definition of the magnitude 
of a vector. 

At this stage it is necessary to introduce the contra- 
variant tensor conjugate to which can conveniently 
be written g*K Then equations (18.1) read 

(16.2) g„g<* = 
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We now define the magnitude B of the covariant 
vector B { by the equation 

(16.8) {B)* = e (B )g i) B i B J , 

where is the indicator of the vector B { , and it is clear 
that £ is an invariant. 

A vector, whose magnitude is unity, is called a unit 
vector. It follows from (15.8) that dx i jds is a unit con- 
travariant vector. If the magnitude of a vector is zero, 
it is called a null vector. The tangent vector to a null 
curve is a null vector. 


§ 17. Associate tensors 

The inner product of the fundamental tensor g (i and 
the contravariant vector A* is the covenant vector g^A*, 
which is said to be associate to A f . We define 


(17.1) A, = g„AK 

Similarly we define 

B‘ = 

and say that the vector B { is associate to the vector B t . 

The relation between a vector and its associate is 
reciprocal, for the vector associate to A { is 

g**A, = g”g ik A k = d* A* = A *. 

This process of association is often referred to as ‘lowering 
the superscript’ or 'raising the s ubscript’ respectively. 
We have 


e {A) g {t A i A* — e (A) g ii g <k A k g f, A l = e (A) g u A k A lt 

which shows that the magnitudes of associate vectors 
are equal. 

The process ofraising and lowering indices canbe perfor¬ 
med on tensors. From the tensor^*, m we can form associ¬ 
ate tensors like or =g ir g l *g m ‘^ w * m . 
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The dot notation is introduced to indicate which indices 
have been raised or lowered. The dots will be omitted 
when there is no possibility of confusion. For example, we 
shall write A** = g ir g is A ra . Note very carefully that 
although g {i and g if are conjugate tensors, the tensors 
A ti and A if are not as a rule conjugate. 

Ex. Show that (A)* = e(A)A f A f . 

§ 18. Angle between two vectors - orthogonality 
The angle between two unit vectors A 1 and B* is defined 
by 

(18.1) cos 0 = g {j A*B* = A,Bi = = A k B k . 

This equation will be found to agree with the usual for¬ 
mula cos 0 as U' + mn' + nn' for the angle between the 
unit vectors (l, m, n) and (V, m', »') in a three dimensional 
Euclidean space when it is referred to rectangular car¬ 
tesian coordinates. We will now show that (18.1) always 
defines a real angle between two real vectors if the metric 
of the Riemannian space is positive-definite. For then the 
magnitude of the vector XA { + pB* is greater than or 
equal to zero for all real values of A and p. That is 

g„(A^< + pB*){XA* + pB>) ^ 0, 

which reduces to 

A® -|- 2A/x cos 0 -(- /t® ^ 0. 

That is 

(A -f p. cos 0) 2 /i a (l — cos 2 0)^0. 

Since this equation is true for all values of A and p, 
it follows that 1 — cos 2 0^0. Thus | cos 0 | ^1 which 
means that 0 is real. If the metric is not positive-definite 
then the angle between two real unit vectors need not 
be real. 

An immediate deduction from (18.1) is that the angle 0 
between two vectors A* and B { , which are not necessarily 
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unit vectors, is given by 


(18.2) 


cos 0 = 


Ve {A) e m g ln A l A”'g n B'B'' 


Two vectors are said to be orthogonal to one another 
if the angle between them is a right angle. From (18.2) 
the necessary and sufficient condition for the ortho¬ 
gonality of the two vectors A { and B* is that 

(18.8) guAW = 0. 

We do not define the angle between two vectors if one 
or both of them happens to be a null vector, but we shall 
take (18.8) as the definition of orthogonality of two null 
vectors. It follows that a null vector is seif-orthogonal. 


Ex. Prove that (1, 0, 0, 0) and (v% 0, 0, V 8 / c ) are unit vectors in 
the V t with the metric (14.2). Show also that the angle between 
these vectors is not real. 


§ 19. Principal directions 
From the symmetric covariant tensor A (j we can 
construct the determinants! equation 

(19.1) \A (i — Xg if |=0, 

which is of degree N in A. When we change to a new 
coordinate system z ( , this equation transforms to 


C^lk — A £lk ) 


dx l dx k 
dx * dx 1 


= 0. 


On applying the multiplication law of determinants, this 
equation can be written 

I Aik — Aglk 


dx 1 
dx { 


= 0. 


The Jacobian determinant 


dx 1 


dx 1 


does not vanish; con¬ 


sequently this equation reduces to 

I An ;—| = o. 
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By comparison with (19.1), we deduce that the roots 
A(/r) of this equation are invariants. (Capital letters are 
used to label the roots. Their enclosure in brackets em¬ 
phasises that they have no tensorial significance, and the 
summation convention is not to apply to them.) 

Now consider the N equations 

(19.2) (^ o -A m g w )Lfo-0, 

where A(jq is a simple root of (19.1). These determine 
the ratios of the N values of We cannot deduce 
immediately from the quotient law that is a contra- 
variant vector, because the tensor A (i — g<y, with 
which it shares inner multiplication, is not arbitrary. 
Instead we change to the coordinate system x { , and equa¬ 
tions (19.2) transform into 

(A lk — X iK) g lk ) faf l \k) = °- 

Inner multiplication by dx^dSS" 1 yields 

(■^Im Elm) L\k) ~ 


These N equations determine the ratios of the N quan- 
. , dx l 

titles fa; L\k) » which are the components of in the 

x* coordinate system. That is, transforms in accor¬ 
dance with equations (5.1) and thus V (K) is a contra- 
variant vector. Let us now choose the ratios of the 
so that it is a unit vector. That is, 


(19>8) gt/L^L^ — cjff), 

where is the indicator of the vector L\k)- 

Similarly we can show that corresponding to each 
simple root A.(m) of equation (19.1) there corresponds the 
unit vector satisfying the equations 

( 19 * 4 ) (A if — ^(U)gn)L\u) — 0 
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and 

(1®’® ) Sn^\u) L\m) ) = £(*o • 

When A ( m) is a multiple root of (19.1), the vector L\ i{) 
is not uniquely determined by equations (19.4) and (19.5). 

Let us choose two simple roots o and A(jq of (19.1). 
Since A(jo # A(m)» the inner multiplication of (19.2) by 
L U and (19.4) by gives us on subtraction, the 
equation 

(19.6) = 0. 

This shows that the two unit vectors L \k\ and L { {U) are 
orthogonaL Thus if all the roots of equation (19.1) are 
simple at any point, the covariant symmetric tensor deter¬ 
mines uniquely N mutually orthogonal unit vectors. The 
directions of these vectors at a point are called the 
principal directions determined by A {j . If the metric 
g ij dx i dx i is positive-definite all the roots of (19.1) are 
real*. That is, the principal directions are real in a space 
with positive-definite metric. 

Now consider the invariant X defined by 

_ A fi UL* 

*~g tm L'L» 

The finite maxima and minima of X are given by = 0, 
that is, by ^ 

A ( ,L<(g tm L'L’»)-g ii L<(A lm L'L'*) - 0 . 

This equation can be written 

M</ — — 0. 

Elimination of the & yields the determinantal equation 

I Aa — tyu 1=0. 

Thus the finite maxima and minima of A are those values 
corresponding to the principal directions determined 
by A {i . 

♦ W. L. Ferrar, Algebra, p. 145. 
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If A u — Agu at a point, then the principal directions 
are indeterminate at that point. If A u = Xg t) at all points 
of a V N , the space is said to be homogeneous with 
respect to the tensor A {j . 

In a Euclidean space of N dimensions referred to 
rectangular cartesian coordinates, the components of the 
fundamental tensor g u form the unit matrix. Hence the 
roots of (19.1) are in this case the latent roots* of the 
matrix A i} and the principal directions are those of the 
latent vectors. 


Ex. Prove that AyL^Lfa = e (K) A (Jr) 
(No summation over (2£))* 


and 


a u l, \k) l \u) 


0 . 


Examples 

1. Show that the angle 0 between the vectors A* and B‘ is given by 

sin*# = < e M) e (B)gt<g/* —g»g u )A k A*B*& 

e(A)e(B)g t>i &it i A k A l BiB t ‘ 

2. If A {) is a skew-symmetric covariant tensor, prove that (A a /Vg, 
A„/Vg, A lt fVg) are the components of a contra variant 
vector. Show also that (VgA**, VgA ”, VgA 1 *) are the com- 
ponents of a covariant vector if A** is a skew-symmetric 
contra variant tensor* 

8. Prove that no relation of the type 

AgijAh* + figiiAjt + vgaAu = 0 

can exist in a Fjy, (IV > 1), where A, /i and v are invariants 
and ^4^ is a symmetric tensor* 

If ^hk ,s a skew-symmetric tensor, show that this equation 
not exist in a Vs for which N > 2. 


* A. C. Aitken, Determinants and matrices, p. 78* 
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§ 20. Chri8toffel symbols 
Although we found in section 5 that dx { /du is always a 
contravariant vector, the exercise at the end of that sec¬ 
tion should convince us that its derivatives (Px { jdu" do not 
form a vector, whose components in any other system are 
the corresponding second derivatives. Again, in section 6 
we proved that the partial derivatives of an invariant 
form the components of a covariant vector; yet we shall 
show in section 22 that the derivatives of a vector do 
not form a tensor whose components in any other system 
are the corresponding derivatives of the transformed 
vector. Our aim is now to build up expressions involving 
the derivatives of a tensor, which are the components of 
a tensor. In order to carry out this programme we must 
first investigate two functions formed from the fun¬ 
damental tensor g {i . These are the Christoffel symbols of 
the first and second kinds defined respectively by 


( 20 . 1 ) 

and 



1 

Ail 

to i 

dgq\ 

dx i 

dx k ) 


(20.2) { ?•} = g' k lih *]. 

Although we shall see that the symbols [ij, k] and 

are not tensors, this notation is introduced in keeping 
with the summation convention which generally applies 
to two indices, one a superscript and the other a subscript. 
Thus all but one of the indices of the Christoffel symbols 
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are regarded as subscripts. The exception is the l in the 
symbol of the second kind, which is treated as a super¬ 
script. 

The definitions show that both symbols are symmetric 
with respect to the indices i and j. Inner multiplication 
of (20.2) by g tm yields 

(20.8) [if, m] = gm {'.). 

It follows immediately from (20.1) that 

(20.4) fa* = [ij, 1c] + [kj, *]. 


We now wish to express the derivatives of g <k in terms 
of the Christoffel symbols, and so we differentiate equa¬ 
tion (16.2) with respect to x l , and obtain 

Inner multiplication by g im gives us 

-1_ aim ait fa* _ q 

dx l ' ® g dx l 

Substituting from (20.4) and (20.2), we finally obtain 
(80.8) {*}-«**{”}• 

We shall now deduce a useful expression for | M. Diffe- 

\ l l) 

rentiate the determinant g = \ g if \ f remembering that 
g lm g is the cofactor of g lm in this determinant, and obtain 


dg 


— aim 0 

UiV' 

From (20.1) and (20.2) and the symmetry of g if we have 


3 g 
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/ * } = JL gtm ftg m , d&m __ JLgfm Oils , 

\iif 2 s \9a> ^ dx* da?*) 2 8 dx* 

Thus 

(2M) {^}=H^ = ^ (1o8 ^ ) - 

Since g is not an invariant, it does not follow that 

{ .*. I is a covariant vector. If g is negative, equation (20.6) 

V) (i) d __ 

should be altered to | „ j = ^ {log V —g}. 

Ex. 1. Calculate the Christoffel symbols corresponding to the 
metrics 

(a) ds* = (da: 1 )* + (a: 1 )*(<&*)* + (a?)* sin* a*(<fc»)». 

(b) da* = (dir 1 )* + Q(a?, x*)(dx*)*, where O is a function of 
a? and a;*. 

Ex. 2. If the metric of a V N is such that g u = 0 for i show 
that 

W UiJ 2g„ "bn* 

21 “ =T7, I «) - A ( lo 8 Vg^}, 


where i $ j and k are not equal, and the summation conven¬ 
tion does not apply. 


§ 21. Transformation law of Christoffel symbols 

The fundamental tensor g it , being covariant, trans¬ 
forms according to the equation 


( 21 . 1 ) 


_ _ dx 1 dx* 

Slm ~W l dff* gii ' 


On differentiating this equation with respect to a?*, we 
have 


dg lm _ dx* dx* dg it dx k d 2 x { dx* dx* d a x* 
dx* ~ dx 1 dff* dx* d& + Wd& d&* gtt + W dx m d£* gii ' 
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We subtract this equation from the sum of the two 
similar equations obtained by cyclic interchange of the 
indices l, m and n, and divide by two. Then with appro¬ 
priate changes of dummy indices we have 


/<n o\ r?- 1 r** n 3®* 3** 3®* . dx* d a X* 

(21.2) [to, »] = [»,. i] ^ ^ + e „ _ _ 


dx l dx m * 


where the bar over the Christoffel symbol indicates that 
it is calculated in the coordinate system x { with respect 
to its fundamental tensor g i{ . Now the transformation law 
of the contravariant fundamental tensor is 


( 21 . 8 ) 


?' = g 


n dx? dx p 
da? dx * 


Inner multiplication of both sides of equations (21.2) by 
the corresponding sides of (21.8) yields on reduction the 
relation 

(21 41 / P \ _ / 8 \ §•£ 3®* , 3®^ 3 8 ®* 

\ij ] da? dx 1 dx m da? dx 1 dx m ' 

Equations (21.2) and (21.4) comprise the transfor¬ 
mation laws of the Christoffel symbols, and clearly in¬ 
dicate that they are not tensors. However, in the very 
special case of linear transformations of coordinates, 
d 2 x*/dx l dx m = 0 and the symbols then transform like ten¬ 
sors. Now inner multiplication of (21.4) by dafjdx- p yields 

(21 51 3 8 ® r _ f p 1 daf (r 1 dx* dx f 

dx l dx m ylm) dxP [ ij J dx 1 dal" 1 

This important equation expresses the second partial 
derivatives of the x r with respect to the x» in terms of 
the first derivatives and Christoffel symbols of the 
second kind. 


Ex. Prove that the transformations of Christoffel symbols form 
a group. 
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§ 22. Covariant differentiation of vectors 
Let ns investigate the tensor character, if any, of the 
partial derivatives of a contravariant vector. We start by 
differentiating the transformation law 

( 22 . 1 ) A * = A( aTi 


with respect to cA, and obtain 

BA k _dA i dS n dx k , Si d*x k 
9®* dt&* dx* dx l 9x < 3x n 9®^ 

The presence of the last term on the right-hand side of 
this equation shows that the partial derivatives dA k jdx i do 
not form a tensor. To obtain a tensor involving the 
partial derivatives we eliminate the partial derivatives 
of the second order by means of (21.5) and this gives us 

dA* _ dA* 7,5? nlT) _| ft 1 9®^ d®* ~| 

W ~ d& dxi W* + A 9a* [_(»») B& \rs) d& 9a?»J* 

In virtue of (22.1) and by appropriate changes of dummy 
indices, this equation reduces to 

BA k , (ft\ j. _ r BA* , /Tl r r 1 0^ dafi 

dx* ^ (ryj [ d& (rnj J dx? dSf * 

Introduce the comma notation 



so that the above equation can be written 

ri 3^ 9®* 

A *>- ds 

Hence from (8,8) it is obvious that A* f is a mixed tensor 
of the second order and it is called the covariant deri¬ 
vative of A k with respect to a?*. 
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In order to set up the corresponding entity for covariant 
vectors, we may conveniently start by differentiating the 
transformation law 

. , T dx* 

( 22 . 8 ) A { = A iWi 

with respect to x l . This yields 

dA t dA } dx l dx i . d 2 x i 
W~d^dS l dx t+ i d&dx 1 * 


Again, by (21.5), we eliminate the partial derivatives of 
the second order. Further we change dummy indices as 
required and we have on substitution from (22.8), that 

dA { fml j _T dAj (r) .‘] dx / dx? 

dx l \t7) m — \ dx n l/nj r jdx i dx l 

The comma notation 


(22.4) 



A r 


is now introduced, and the above equation can be written 

j _ dx f 9®" 

A t ,i = 


which shows that A Jt „ is a covariant tensor of the second 
order, called the covariant derivative of A f with res¬ 
pect to x n . 

In a Euclidean space of N dimensions, referred to 
rectangular cartesian coordinates, the components of the 
fundamental tensor g i} are zero except g u = g 22 = . . . = 
= g NN = 1. Thus all the Christoffel symbols are zero and 
so covariant differentiation reduces to the familiar partial 
differentiation. It is well to observe that the Christoffel 
symbols do nof all vanish in a Euclidean space referred, 
for example, to spherical polar coordinates. 

We can construct the invariant Aj by contraction. 
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Applying (20.6) we obtain 

4-SMJ}'- 


dAi + 


dx* 


That is 


< 22 - 5 > A ‘< - 

This invariant is called the divergence of the contra- 
variant vector A* and is often denoted by div A { . The 
divergence of a covariant vector A t is defined by 

(22.6) diV A f = g* k A itk . 

The partial derivatives of an invariant form the com¬ 
ponents of a covariant vector. We extend the definition 
of covariant differentiation to invariants by calling the 
familiar partial derivative the covariant derivative. That 
is, from the invariant I we form by definition 



Since I { is a covariant vector, we find from (22.4) that 
its covariant derivative with respect to cA is given by 


(22.7) 


,,, _ jy_fri ar 

' ,t ' ,i dx i dx i \ ij j dx T 


Thus (I *) * = That is, the covariant differen¬ 

tiation of invariants’is commutative. In section 81, we 
shall see that the covariant differentiation of vectors is 
not in general commutative. 

We can form the divergence of the covariant vector 
1 1 % this is called the Laplacian of I and written P 8 !. Then 


r*i = g ,k (i.,U = s'* 



Ex. 1. Show that A itn — A n , f 


*bA, 
IV x n 


lx* ' 
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Ex. 2. Prove that div A f = -iL — = div A*. 

Ex. 8. Evaluate div and V*I in (a) cylindrical polar, (b) spherical 
polar coordinates. 


§ 28. Covariant differentiation of tensors 
In the last section we constructed tensors containing 
the partial derivatives of vectors. Can we now extend the 
process of covariant differentiation to tensors? For our 
typical tensor we shall choose A}. There will be no loss 
in generality involved as this tensor has one contra* 
variant and one covariant index. Inner multiplication of 
its transformation law (8.8) by daP/dx* yields 


(28.1) 


11 

A -w 


= Af 


dx l 
dx* ’ 


We differentiate with respect to x k , then eliminate the 
second order partial derivatives by means of (21.5) and 
obtain 


dAj dx m ii[\P\ d-g” 1 jfar da ?"j 

dx k dx* f [_[ ik | dxP \rs J dx* dx k J 

__ dAf fa:* dx l Rpl 9®*_ j l\dx T da? 1 

dx* dx k dx* 1 |_\ jk J dx p \rs] dx* dx k J* 

Applying (28.1) to this and changing appropriate dummy 
indices, this equation becomes 



Introduce the comma notation 


dx* dx* 
dx k dx* ’ 


(28.2) 


A m — 
J1 l,t = 


dAf 

dx* 


+ 



Af, 
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and the above equation takes the form 

n _ Am 3»* dx l 

A ‘*W ~ A ''*d&d&' 

On inner multiplication by dx r Jdai m we have 

jr _ .m d£r dx ' dx * 

i,k l,t dttf" dx l dx k * 

Hence A™ t is a tensor of the third order of the type 
indicated by its indices. This tensor is called the co¬ 
variant derivative of A f with respect to x*. 

An examination of equation (28.2) shows that the 
covariant derivative of Af contains three terms. They 
are (1) the partial derivative, (2) a term with positive 
sign similar to that contained in the covariant derivative 
of a contravariant vector and (8) a term with negative 
sign similar to that contained in the covariant derivative 
of a covariant vector. This suggests that the expression 




(28.8) 


daf 1 


"!p 


+ Z (M A u i" u 
^txlfcnj r i" r « 

Ir 

J r t** r 0-l * f j 


P 

— E 
0=1 


r fi n I 


is a tensor, called the covariant derivative of A“ l "“‘ 
with respect to a". The proof which we worked out for 
the mixed tensor of the second order applies also to (28.8). 
We will omit its tedious details and supply an alternative 
proof in section 80. 

On referring to equations (20.8), (20.4), (20.5) and 
(28.8) we deduce that 


(28-*) g«.» — {,4} {,*} g « ~ °’ 
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and 

<*»•«> «.-J3+{i} 

The covanant derivatives of covariant derivatives are 
again tensors. We indicate these covariant derivatives 
of the second order by adding another index without a 
comma. For example A itfk is the covariant derivative of 
A tli with respect to x k . 

Ex. 1. If A u = B t>/ — B,,„ prove that A„, k +A n , t + A k{ ,, = 0. 
Ex. 2. By means of (28.5), show that div A* = div A { . 

Ex. 8. If A i)h is a skew-symmetric tensor, show that 
1 d 

—y r— : WgA***) is a tensor. 

■vg 

§ 24. Laws of covariant differentiation 
Covariant derivatives obey the following laws:- 

(1) the covariant derivative of the sum (or difference) 
of two tensors is the sum (or difference) of their covariant 
derivatives. This law is an immediate deduction from 
(28.8). 

(2) the covariant derivative of an outer (or inner) 
product of two tensors is equal to the sum of the two 
terms obtained by outer (or inner) multiplication of each 
tensor with the covariant derivative of the other tensor. 
Consider as an illustration 

= + A" B l m . 

This type of proof is quite general and will apply to any 
case of outer multiplication of two tensors. (Another 
proof is provided in section 28). Contraction of l and j 
gives us 

+ A„0„, 

which makes it clear that the rule also applies to inner 
products. 
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(8) the tensors g ijt g iS and 8j are constants with respect 
to covariant differentiation. This is merely another way 
of stating equations (28.4), (28.5) and (28.6). For example 

(S {i A,). m = CAi + = g u A iUn . 

Ex, If I and J are invariants, show that 

div = JV a I + gf'I.iJ.r 


§ 25. Intrinsic derivatives 
Consider the tensor A“ 1 '‘“‘ whose components are 
functions of t along a curve x* = x { (t). The intrinsic 
derivative is defined by 


(25.1) 


81 ri " r «” 


dx k 


Accordingly, the intrinsic derivative is a tensor of the 
same order and type as the original tensor. 
Corresponding to the invariant I, we have 

81 _ j dr*_a/dr*_dZ 
8t '* dt ~ 9®* dt — dt * 


That is, the intrinsic derivative of an invariant coincides 
with its total derivative. 

Intrinsic derivatives of higher order are easily defined. 
For example, 

8*A< _ 8 (8A<\ ( { dx k \ 

8t? ~ 8t\8t) ~ dt J., dt ’ 

In general, intrinsic differentiation is not commutative. 

From (22.2), (22.4), (28.4), (28.5) and (28.6) we cal¬ 
culate that 


(25.2) 

(25.8) 


Ad* _ dd* m 

dt ~ dt + \r?j 
8A k dA k (r) 

~8t = ~dt \jbyj 


A dt ’ 
. dx s 

d 'W 
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tea 

dt 


dg” _ 36} _ 


6t 


dt 


From the definition of intrinsic derivatives, it follows 

that they obey the same three laws which apply to 

covariant derivatives, 

(5 /daj<\ d*x { l i \ dx f dx* 

Ex, Show that — (— 1 = —— + {..} — — • 
dt\dt) dt* dt dt 


p. 28 Ex. 1 
kind are 

(a) 


Solutions 

The only non-zero Christoffel symbols of the second 


1 

22 

2 

88 

i* 

\22 




{m} - — *** to,, “*- {, 2 ) - 1 ^‘- 

’*'■ (is) ■ ,,x '' (m) ■ 


= — sin x* cos - 


v / 

_1 IG 

~~ T jj® 3 ’ { 12 } ” 2 G lx' 


(b) 

p. 88 Ex. 8. 

, . .. .. IA' IA* IA» 1 

(a) div A* = —- + r-T + ITT + 1 A * 
da 1 bx* da? 3 x l 


1_IG 
2 G 7>x 3 


V*I = 


bA x 

bx l 

b*I 


1 ^ 4 , ^ 4 , 1 

(x')' lx' + lx* + a 3 " 


((5a: 1 ) 2 (a; 1 ) 2 (baj a ) a 


1*1 * 3 I 


_+IH 

(t)® 3 ) 3 x' lx' 


. „ IA* d.4 3 2 .... 

(b) div A' = — + —— + — + A 1 + cot x 3 A a , 


bx l 

lA t 


bx* 

1 


bx* 


_ + _m- + _;_ 

<)a? 1 (a; 1 ) 8 bx 3 (a; 1 ) 3 sin 3 a? 8 bar* 


IA, 2 t»t* 3 

sST + 3 A * + — A « 


v*i = 


x 1 - (a;!) 3 

b*I 2 bl . cot x* bJ 


b*I 1 b*I 1 d 3 / 2 bI cot a; 3 d/ 

(bx 1 )* + (a! 1 ) 3 (i>® 3 ) 3 (a: 1 ) 3 sin 3 ® 3 (<5x 3 ) 3 + ®> ii® 3 ^ (a 1 ) 3 da: 2 


The corresponding results in Rutherford’s Vector Methods, pp. 
72—78, differ slightly because they are expressed in terms of the 
physical components of the above vector (see section 62). 
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§ 26. Geodesics 

In Euclidean three-dimensional space, a straight line 
is the path of shortest distance between two points, and 
it is our aim to generalise this fundamental concept to 
Riemannian spaces. Let C be the curve X 1 — x*(t) with 
parameter t joining two fixed points P 0 and P v whose 
parameters are t Q and ^ respectively. Then the distance 
8 along the curve between P 0 and P l is given by 


/ou i \ f‘ l 1 / dx* dad 

<SM) S= J,. m dL 

Consider all the curves passing through the two fixed 
points P 0 and P v Any of these curves, for which the 
distance P 0 P X measured along the curve is stationary, 
is called a geodesic. We could obtain the differential 
equations of the geodesics by applying Euler’s equations, 
a well-known result in the calculus of variations, to 
equations (26.1). However we shall find it more instruc¬ 
tive to appeal to first principles. 

Let us choose a small arbitrary vector dx { varying con¬ 
tinuously along C. Then the equations = x l -f da* 
define a neighbouring curve C to C. Further let us impose 
the conditions that dx i = 0 at P 0 and P v This means 
that the curve C always joins P 0 to P v The distance s 
from P 0 to P 1 along C is given by 



dx* dxi 

es " w *■ 
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where the g tJ (x) are now functions of x l . We have 


...dx* dx { _[ dg„ x Alda* 
gii ^ dt dt~ ( gw + ' 


= Su 


dx* 
dx i dx i 


dt 


d(dx*)\(da* d{6x*)\ 

dt ' dt )\dt + dt ) 
da* d(6x*) dg ti a _ b da*da^ 
dt +2gi ‘ dt dt + ** dt dt ’ 


where terms of higher order than the first are neglected. 
Thus 


t 

t 


. dS* dx } 


dx i dec* 

egt, HtHt 


1 + 


dx { d{6x s ) 1 dg if , da* da f ~\ 

e "-5-~5~ + ~5"S 

dx i dot* 
gii ~dtlt 


Consequently the variation in length 6s from the curve 
C to the curve C is given by 

5s = 8 — s 

p*i da* d(6x?) 


gii ~dt 


dt 


JL — — 

^ 2 da k dt dt 


4 <„ 


1 / da* da 1 
r egii ~dt ~dt 


dt. 


We now simplify this equation by choosing the arc- 
distance s along C as parameter and obtain 

r*i r da* d(6x 0 1 dg { , a . da { daH ds, 

61 '-J„ ['“a— +1 

where s 0 and s x are the values of s corresponding to the 
points P 0 and P t respectively. Integration by parts yields 

The integrated portion vanishes since 6x j is zero at P 0 
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and P v Also we have 
d / d®*\ _ d?x* dgij dx* das* 

ds\ 8ii ls)~ S{i ~dF + dx*Ha~ds 

cPx* , 1 dg if dx* dx* , 1 dgki dx* dx* 
~ 8ii IF + ~2 d&~da 1* + t ~dx* Is ~ds 


It follows that 


f*t r d?x* ,,, dx* dx*~\ 

(28.2) i, -£ M [g„ W + W. >] * -5-] *• 

The variations <5®* are arbitrary, thus the necessary and 
sufficient conditions that the curve C be a geodesic are 

d*x* dx* dx k 

(26-8) = 0. 


Inner multiplication by g il yields the contravariant form 


(26.4) 


d_ (—\ = — f M — Act 

as \ds/ “ ds» + \iJfcj ds ds 


Either set of equations (26.8) or (26.4) are the differen¬ 
tial equations of a geodesic. They constitute N differen¬ 
tial equations of the second order. The theory of differen¬ 
tial equations states that a solution x* — x*(s) is deter¬ 
mined uniquely if the initial values of x* and dx*fds are 
given at any point. Geometrically this means that there 
is a unique geodesic with given direction at any point 
of the space. We defined the geodesic in terms of the 
curve passing through two points, but this geodesic may 
not be unique unless the two points are sufficiently close 
to one another. The problem of uniqueness now involves 
topological properties of the space V w . For example, there 
is a unique geodesic passing through two points on a 
sphere, except when the two points are at the ends of 
a diameter. In this latter case, all great circles passing 
through the two points are geodesics. 

For a Euclidean space, referred to rectangular cartesian 
coordinates, the Christoffel symbols are zero. Hence the 
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geodesics are given by d 2 x l /ds 2 = 0, whose solution is 
x l = A l s + where A 1 and B l are constant vectors. 
That is, the geodesics are straight lines. 


§ 27. Null-geodesics 
Equations (15.8) state that 


(27.1) 


dx { dx* 

g " ds ~ds = 6 


along any portion of a curve which is not null. On diffe¬ 
rentiation we obtain 


d 

( dx i daA 

d 

/ dx i dx<\ 

dx { 6 (dx f \ 

ds 

[ gii H'fo) 

“ds 1 

ds ds) 

— gii ds ds\ds) 


It follows from (26.4) that the invariant ^ 

is zero at all points on a geodesic. Thus the indicator e 
cannot change abruptly along a geodesic, and so if the 
tangent vector is not null at any one point, it cannot 
be null at any other point on the geodesic. On the other 
hand, if the initial direction is null, then the curve is 
null and it is of course impossible to introduce the arc- 
distance as parameter. Instead we now say that a null 
curve x i — **(/) which is a solution of the equations 



(27.2) 


d 2 x l t n dx* dx k 
dt 2 + (ifrj dt dt 


is a null-geodesic. 

The null-geodesics in the V x with line-element (14.2) 
satisfy the equations = 0. Therefore the null curves 


given by (15.2) do not satisfy these equations unless r, 
6 and ip are constants. It follows that a null curve is not 
necessarily a null-geodesic. 


§ 28. Geodesic coordinates 
We shall now show that it is always possible to choose 
the coordinate system so that all the Christoffel symbols 
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are zero at a particular point. Consider a general coor¬ 
dinate system x ( , whose values at a particular point P 0 
are x* 6) , and introduce a new coordinate system x* by 
the equations 

(28.1) 0 = + i {j n ) B («”—*5,). 

The index (0) attached to any entity denotes its value 
at the point P 0 . The brackets serve to emphasise that 
this index has no tensorial significance and that the 
summation convention does not apply to it. Differentia¬ 
tion with respect to a* yields 


(28.2) 


Hence 


d& 

dx i 


©„-*> 


= <3j+ (.*’) 


(®» 


U (0) 


)• 


Accordingly the Jacobian deter- 


minant 



is not zero, which shows us that the 


transformation (28.1) is permissible in the neighbour¬ 
hood of P 0 . On inner multiplication of (28.2) with 
dx i /dx^, we obtain 

dx* i i \ . n dx* 

k ~ dx k |jn)(o) 3,(01 fa*' 

We differentiate this with respect to x 4 and obtain 

•H&+W.55+KL*-*-- ,w 

Thus at P 0 

(&)(o, = (dl^lco) =_ {/n} (0) = - {A} 

We now substitute in (21.4) and derive 

W« " {*?}«» 8 ‘ 6 ^ — ^ (lm) (0) ' 

That is, ,— 


d^dx 4 ’ 
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Hence a particular system of coordinates, called geodesic 
coordinates, can always be chosen so that the Chris- 
toffel symbols are zero at any assigned point called the 
pole. The transformation (28.1) is not the only method 
of obtaining geodesic coordinates. However, we have 
shown that one such system exists and in this particular 
system — 0, which means that the pole is also the 
origin of coordinates. 

The important property possessed by a geodesic coor¬ 
dinate system is this: the covariant derivatives reduce to 
the corresponding partial derivatives at the pole because 
at this point all the Christoffel symbols are zero. We men¬ 
tioned earlier in section 9 the fact that if a tensor is zero 
in one coordinate system, it is zero in every coordinate 
system, since the transformation law of tensors is linear. 
The setting up of a tensor equation often involves heavy 
algebraic manipulations. Generally, the volume of work 
is reduced by first proving the equation with respect to 
a geodesic coordinate system at its pole. It follows that 
the equation is true for all coordinate systems at this 
point. Then if this point is general, the equation is true 
at all points of the V N . 

We shall illustrate this method by proving the multi¬ 
plication law of section 24 which is satisfied by covariant 
derivatives. Consider the tensor 


Let us choose a geodesic coordinate system with pole at 
then the covariant derivatives reduce to the familiar 
partial derivatives at P 0 . Since partial derivatives satisfy 
the multiplication law 


d_ 

du 


<"> = In*’+ ’’§;■ 


the above tensor is zero at P 0 in the geodesic coordinate 
system. Thus it is zero at P 0 in every coordinate system. 
But P 0 is a general point, therefore the above tensor is 
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zero at all points in the V#. This establishes the truth 
of the multiplication law. 

Ex. Show that at the pole P 0 of a geodesic coordinate system 

a » /0 

* i)** (ij j * 


§ 29. Parallelism 

An important property of parallelism in a Euclidean 
space, which is referred to rectangular cartesian coor¬ 
dinates, is this: a parallel field of vectors A t is obtained 
throughout a Euclidean space if the components A ( are 
constants. We can express this analytically either in the 
form dAf/dt = 0 or by dA i jdx i = 0. Since the Christoffel 
symbols are zero, we can write these equations equi¬ 
valently in the tensorial forms dAJdt = 0 or A {ti — 0 
respectively. This suggests two ways of generalising the 
concept of parallelism to a Riemannian space. We shall 
see in section 81 that the partial differential equations 
A tij = 0 are in general not consistent. Thus the second 
suggested generalisation is not a profitable one. Also the 
intrinsic derivative dAJdt is only defined along a curve, 
hence, using the first suggestion, we can only define 
parallelism along a curve. Formally, the vectors A { 
constitute a field of parallel vectors along the curve 
x { = x*(t) if A t is a solution of the differential equations 



These equations form a set of N differential equations 
of the first order, and consequently if the vector A { is 
given at any one point on the curve, it is uniquely deter¬ 
mined at all other points of the curve. We may also say 
that a field of parallel vectors is obtained from a given 
vector by parallel propagation along the curve. Since 
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T-S 


<5* 


we can write the condition for parallelism along a curve 
in the contravariant form 


(29.2) 


dA* _ dA' 
St dt 



= 0 . 


We see from equations (26.4) that the unit tangent 
vectors form a field of parallel vectors along a geodesic. 

The magnitude A of the vector A i is given by 
{A) z — e (A ) g ij A i A i . On differentiation we have 


.dA d . d f\A* 

2A ~dt~dt ( e W&i A Ai ) = -fo (e ( A)g ii A i A^ = 2e {A) g ( , — AK 


This equation becomes A(dA/dt) = 0 if A 1 forms a field of 
parallel vectors. We deduce that the magnitude of all 
vectors of a field of parallel vectors is constant. 

The angle 6 between the two vectors A* and B* is given 
by AB cos 6 = g {j A i B i . We obtain on differentiation 

AB± (cos 6) + B+A i?)cos 6- Sll 

When both A { and B { form fields of parallel vectors, this 
equation reduces to d(cos d)fdt = 0, provided that 
neither of the vectors A i and B* is null. That is, the 
angle between two non-null vectors remains constant 
whilst both undergo parallel propagation along the same 
curve. 

The vector obtained at Q by parallel propagation 
from P depends on the curve joining P to Q. Hence 
parallel propagation around a closed curve does not neces¬ 
sarily lead back to the initial vector. As an example con¬ 
sider the V 2 formed by the surface of a unit sphere. On 
choosing spherical polar coordinates x l = 0, a? 2 = y the 
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metric of the V i is given by ds 2 = dd 2 + sin 8 6dy> 2 and a 
brief calculation shows that the only non-vanishing Chris- 


toffel symbols of the second kind are I = — sin 6 cos 6 

(2 1 l 22 / 

and | | — cot 6. Let us select the small circle 6 = a for 

the purpose of parallel propagation of the vector A 1 . 
Along this circle dQ/dt — 0, and consequently the equa¬ 
tions (29.2) reduce to 


dd 1 . Ja dA 8 

—-cos a sin ccA 2 = 0 -, -=- 1 - cot a A 1 = 0 

dip dip 

whose solution can readily be obtained in the form 


A 1 = sin a [c sin (ip cos a) + d cos (ip cos a)]; 

A 2 = c cos (ip cos a) — d sin (ip cos a), 

where c and d are constants. Suppose we choose the vector 
A { to be (1, 0) at the point defined by ip = 0. Then by 
substitution we have c = 0 and d = cosec a. Thus the 
vector A* is uniquely determined by the components 
(cos [y cos a], — sin \ip cos a]/sin a). The result of parallel 
propagation around the small circle then yields the vector 
(cos [2jr cos a],—sin \2n cos a] /sin a) which differs from the 
original vector (1, 0). The case of the great circle a = ti/2 
is exceptional, as along it parallel propagation does lead 
back to the original vector. 


Ex. In the V t with metric da 1 — du 4 + SS Xdudo + do 1 , where A 
is a function of u and o, show that the tangent vectors to the 
curves ti = constant form a field of parallel vectors along the 
curves o = constant. 


§ 80. Covariant derivative 
By means of the concept of parallelism we are now 
able to supply a proof that the right-hand side of equation 
(28.8) constitutes a tensor. Consider a curve C determined 
by the equations x‘ = x*(t). Choose p arbitrary contra- 
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variant vector-fields Z‘ a) , X{ 2) , . . Xj p) each parallel 
along the curve C and s arbitrary covanant vector-fields 
which are also parallel along the curve 

C. Then 


(80.1) 


. M 

dt ^ \jk 


IF 


= 0 (/? = 1, 2,... p), 


(80.2) 


(«)< 

dt 


= 0 (a «■ 1,2,... a). 


Now let us set up the invariant 

1 — A ri:'.r* x \i) X'd) • • X \p) ^(l)^ ^(a)Uj • • 

The derivative of this invariant with respect to t is also 
an invariant, and by application of (80.1) and (80.2) and 
appropriate changes of dummy indices we find 


dt “ ■* 8 ) * * ^\v) • • Y W u, 


dA? 

_ r i • • r » 

dt 


a-1 r ‘" r<> 


-XA u i" u * ( 1 

pmi r i" r s-i ir B+i - T s\rpn 


\fcraj 


We deduce from the quotient law that the expression m 
brackets on the right-hand side of this equation is a 
tensor, which we call the intrinsic derivative, and denote 
by 8A*iy*'[dt. It follows immediately that 


dAp"? 

r * • • f ' 


“• daf 1 * 

^ = ~ » + ZAp”?*- !*“«« 

dt |_ daF ^ r »" r p 


x r ‘" r * \ftnj 


!• HI HI 



48 TENSOR CALCULUS § 30 

The quotient law again shows us that the expression in 
square brackets is a tensor, which we denote by 
and which we call the covariant derivative. This completes 
the required proof. 



CHAPTER V 


CURVATURE TENSOR 

§ 31. Riemann-Christoffel tensor 

We will now investigate the commutative problem with 
respect to covariant differentiation. Let us begin with the 
covariant derivative of an arbitrary covariant vector Aj 



A further covariant differentiation yields 



(i\dA t iinmai, 
dx n dx p \jnj dx p 1 dx p \jnj \jpj dx n 



We interchange the indices n and p and subtract. After 
changing several d umm y indices we have 

Af,np Aj tPn = 

[» {)!>} 

Since A t is an arbitrary vector, it follows from the quotient 
law that the expression in square brackets is a mixed 
tensor of the fourth order, of contravariant order one and 
covariant order three. Using the notation 


dx p \jnf \jp\ \ps| \jn\ 1 ’ 
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(81.1) B\ 


= _d_ 

inp ~daf' 


1 1 \ d (Mil Ml 4 ! H H s 1 
\jp\ dx v \jn\ \n«j\/pj \ps( ( 771 / ’ 


we observe that R\ inp is a tensor of the fourth order, 
called the Riemann-Christoffel tensor. It is formed 
exclusively from the fundamental tensor and its 
derivatives up to and including the second order. This 
tensor does not depend on the choice of the vector A { . 
We can now write 


( 81 . 2 ) Aj t „„ Aj if> „ — Rj np A t . 

It is clear from this equation, that the necessary and 
sufficient conditions that the covariant differentiation 
of all vectors be commutative, is that the Riemann- 
Christoffel tensor be identically zero. We stated in section 
29 that the equations A t i = 0 are not as a rule consistent. 
In fact, equations (81.2) show us that a necessary, but 
not sufficient condition, for the consistency of A t i = 0 
is that Rj np A t = 0, a set of equations which is not 
generally satisfied. 

Referring to the definition (81.1), we observe that 

(81.8) R] lnp = - R] Svn . 

That is, R\j np is skew-symmetric with respect to the 
indices n and p. 

Ex. 1. Prove that 

Ex. 2. Prove that B l lnp = 0. 


§ 82. Curvature tensor 

We now introduce the covariant curvature tensor 
defined by 

(82.1) Rrinp^grl#,*,- 

On substituting from (81.1) into (82.1) we obtain 
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+ & C,| -!.!,! !J' *» !/„} 

which reduces on application of (20.8) and (20.4) to 
R rin, = ^ &P> r]- 2-. Ijn, r] + { /J [rp, l] - {7 } [m, I], 


We can reduce this further by means of (20.1) and (20.2) 
and finally obtain the important formula 


(82.2) 


(^grp 

, 3‘S,. 


3‘e,. \ 

ydaJdx" 

~ dx r dx p 

dx i dx p 

dx r daf*J 


+ g u ({jn, s] [rp, t] — \jp, s] [m, /]). 


From this result we immediately deduce the relations 


(82.8) 

and 


■®r/nj> = Rjrnp » 

= ^ripn * 


(82.4) R r j nj) + Rfnpj + R rvln = 0. 

We are now faced by this problem: how many distinct 
arithmetical non-vanishing components does the tensor 
R r j„ v generally possess? Referring to (82.8) we see that a 
component is zero if either r = j orn = p. Thus the com¬ 
ponents apart from sign conform to the three types 
Rfirv an d R r jnp w h er e r, j, n and p are distinct from 
one another. There are as many components of the type 
R rirf as there are ways of combining r and j. That is, 
^N(N — 1). There are as many components of the type 
R r irp 88 there are ways of combining j and p after selec¬ 
tion of r. That is, &N(N —1 )(N — 2). The number of 
combinations of r,j,n and p is YtN(N—l){N—2){N—8). 
But R rtnv is determined, except for sign, when r is paired 
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with either j, nor p. Thus there are $N(N—l)(N—2)(N—8) 
components of the type J? r j np . Equations (82.8) enable 
us to rewrite (82.4) in the form 

B imp "(■ B ipm "t" B inpr ~ ®‘ 

And so with a given set of four indices, only one in¬ 
dependent equation of the type (82.4) exists. Further, 
if the indices j, r, n and p are not distinct, (82.4) 
reduces to one of the equations (82.8). Hence the 
number of independent equations of the type (82.4) is 
-^N(N — 1 )(JV — 2 )(N — 8). Therefore the number of 
distinct components of R r)nv is 

%N(N—1)+$N(N—1) (N —2)+ (N —I) (N —2) (N —8) 
— &N(N — 1)(N — 2 )(N — 8 ) = — 1 ). 

In particular we see that the curvature tensor of a F a 
has only one distinct non-vanishing component. 

Ex. 1. Prove the relation (82.4) by setting up a geodesic coor- 
dinate system. 

^ 2(2 

Ex. 2. Prove that B iaia = — G —— for the V a whose line- 

ou a 

element is ds i = du* + G a dv a , where G is a function of u and r. 


§ 88. Ricci tensor - Curvature invariant 

At first sight there appear to be three different ways of 
contracting the Riemann-Christoffel tensor R] inp - We have 
B .*i«p = g UB sinp = °> because R tlnp is skew-symmetric 
in a and l. We see from (81.8) that R l jn , = — R\ ftn . 
Hence we need only consider the contraction, called the 
Ricci tensor, defined by 

(88.1) R in = R\ m = g u R aM - 

On contraction of l and p in (81.1) and on substitution 
from (20.6) we find that 
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R,n = 


( 88 . 2 ) 


dx 1 dx n dx l (/» 


- (A) i f og 


from which it is clear that R in is symmetric. (If g is 
negative, we must replace log y/g by log V—!•) 

The curvature invariant is defined by 

(88.8) R = gin Rjn . 

A space for which R i} = Ig if at all points, where 1 is 
an invariant, is called an Einstein space. Inner multi¬ 
plication by g fi shows that R = NI. Thus for an Einstein 
space 

/ n n . V n I *. 


(88.4) 


— jy RSa- 


Ex. For a V t , prove that gR„ = —g tt R mt and gR = — 2R l31l . 
Hence deduce that every V t is an Einstein space. 

§ 84. Bianchi's identity 

Let us choose a system of geodesic coordinates. On 
referring to (81.1) we have by covariant differentiation 
that 

B ' snvr = h {R '' n * ] = Wd& {jp} - dJfo ( ? i) 

at the pole. Cyclic interchange of n, p and r gives us two 
other equations. We obtain by addition 

< 84 ;1) «!(.*. + + K,.., = «• 

This is a tensor equation, true at the pole of a geodesic 
system of coordinates. Thus it also holds for every coor¬ 
dinate system at that pole. Further any point can be 
chosen as the pole of a geodesic coordinate system. 
Therefore equation (84.1) is true at all points of space. 
Inner multiplication by g lm yields the Bianchi identity 

(84,2) Rfninp t r + = 0. 


Mn/pr, n 


h mirn,p 
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The Einstein tensor is defined by 
(84.8) = — 

The inner multiplication of (84.2) by g mt> g fn and the 
application of (88.1), (88.8), and (82.8) give us the 
equation 

which may be written 

(84.4) R r = 2g»R frtn . 

Differentiating (84.8) covariantly, we obtain 

& )ti = 

That is 

(84.5) & Ui - 0. 

This equation is important in the theory of Relativity. 

§ 85. Riemannian curvature 
From any two vectors A i and B l at a point of a F#, 
we can construct the invariant R rinv A r A n B i B 9 . Let us 
consider what happens if we replace the vectors A* and 
B l by the two linear combinations 

X* = U { + fiB { , Y* = pA* + rB<, 

where A, p, p and t are invariants. A straightforward 
calculation, with the aid of (82.8), shows that 

R rjnt> X'X»Y>Y> = (At — ppY Rfinv A T A n B s B*. 

Thus the expression R rjn9 A r A n B i B v ) which is an in¬ 
variant with respect to coordinate transformations, is 
almost an invariant under linear transformations of vec¬ 
tors. In order to obtain an expression which is also in¬ 
variant under linear transformations of vectors let us 
evaluate 
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= (Mn+f* B »)(M- n +flB”)(pA f ,+TB t ,){pA’>+TB») 

-(XA t ,+{iB v )(pA>+TB*)(XA j +[iB i )(pAi+TB*) 

= (eAX a A a +e B p a B a +2X/i cos 0 AB)(eA p a A 2 +e B x a B z +2pr cos 6AB) 
— (e A XpA a +e B pxB 2 + [Xx+pp] cos 6 AB ) a 
= (Xx—pp) a (e A e B — cos 8 6)A z B a 

= (fo—f*p) i (g rn gjp—grpgjn)4 r A n B i B 9 , 

where 6 is the angle between the vectors A { and B l . It 
follows that 

(85.1) K = - Rrin t , ArAnBi B J> - 

(gmg iv ~ grvgin)A r A n B j B v 

is an invariant which is unaltered at a point, when the 
two vectors determining it are replaced by any linear 
combination. This invariant is called the Riemannian 
curvature of the space V jy associated with the vectors 
A i and 2?*. Note that the denominator of K is unity if 
the vectors A { and J3* are orthogonal unit vectors. 

At any point of a two-dimensional space there exist 
only two independent vectors. Hence the Riemannian 
curvature of a V 2 is uniquely determined at each point. 

Its value is easily found by choosing the two vectors 
whose components are (1, 0) and (0, 1) respectively. Then 

(85.2) K = -^212-_ . 

011022 §12* g 

§ 86. Flat space 

We say that a space is flat if K = 0 at every point of 
it. From (85.1) the necessary and sufficient condition is 

R rinv A r A n BiB> = 0 

for all vectors A* and B*. In virtue of equations (82.8) 
it follows that 

B rjnp "1“ B nin "t” B nvri "t* B rvni ~ 
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That is, 

"t" ^fpni 

which can be written 

^r*np = RrvSn' 

Interchanging j, n and p cyclically we obtain 

■®ro vi ~ ®finp* 

Thus we have 

*®rinp ~ *®rpin 

Substitution in (82.4) immediately yields 

Rfjnt> “ 

Conversely, if = 0, then it is clear that K = 0. 
Therefore, the necessary and sufficient condition that a 
space V N be flat is that the Riemann-Christoffel tensor 
be identically zero. 

In a flat space B] fnp = 0, hence we deduce from (81.2) 
that in a flat space the equations i = 0 are con¬ 
sistent. Inner multiplication by dx*jdt yields dAJdt = 0. 
Thus in a flat space the property of parallelism is in¬ 
dependent of the choice of a curve. We may therefore 
say that parallelism is an absolute property of a flat 
space. 

A familiar example of a flat space is the Euclidean 
plane for which the metric is ds 2 = dx 2 + dy 2 in rec¬ 
tangular cartesian and ds 2 = dr 2 + T 2 d6 2 in polar coor¬ 
dinates. 

Ex. If the metric of a two-dimensional flat space is /(r)[(<£r l ) a + 
(da i 8 ) 8 ], where (r) a = (as 1 ) 1 + (a*) 9 , show that f(r ) = c(r)*, 
where e and k are constants. 

§ 87. Space of constant curvature 
Let us now investigate spaces in which the Riemannian 
curvature at every point does not depend on the choice 
of the associated vectors A i and B*. From (85.1) the 
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necessary and sufficient condition is that 

{K( Srngjp - g rp g in ) - R Tjnp }A r A n B i B p = 0 

for all vectors A* and B*. A similar calculation to that 
of the previous section will show that this condition 
reduces to 

rinp = ^(drnSjji §rj>§fn)» 

where K is now a function of the coordinates x*. 
Covariant differentiation gives us 

SrpSin)’ 

We substitute this result in Bianchi’s identity (84.2) and 
obtain 


^,r(SmnSjp SmpSin) "f* ^ ,n(&mp&ir SmrSfp) 

“I” ^,p(SmrSin SmnSfr ) = 

Inner multiplication by g mn g ip yields 

(N —1)(2V— 2)K r = 0. 

Hence if N > 2, it follows that K is constant. Thus we 
have proved Schur’s theorem; *if at each point of a space 
Vn> (N > 2 ), the Riemannian curvature is a function of 
the coordinates only, then it is constant throughout the 
V N \ Such a Vif is called a space of constant curvature. 

The metric of the F a formed by the surface of a sphere 
of radius a is ds z — a a (dfl a 4* sin 3 ddy> 2 ) in spherical 
polar coordinates. The reader is asked to verify that 
R 1212 = a 3 sin 3 6 and that it follows from (85.2) that the 
surface of a sphere is a surface of constant curvature 1 fa 2 . 

Ex. 1. Show that a space of constant curvature is an Einstein 
space. 

Ex. 2. In a Euclidean V lt prove that the hypersphere 
as 1 = e sin 0 sin <p sin y>, x* = c sin 0 sin <p cos y> 

** = c sin 0 cos <p, x* = c cos 0 

is a F, of constant curvature 1/c*. 



CHAPTER VI 


EUCLIDEAN THREE-DIMENSIONAL 
DIFFERENTIAL GEOMETRY 

Euclidean Geometry investigates the properties of 
/figures which are invariant with respect to translations 
and rotations in space. It may be subdivided into Al¬ 
gebraic and Differential Geometry. The former studies by 
algebraic methods the theory applicable to entire con¬ 
figurations such as the class or degree of a curve. The 
latter discusses by means of the Calculus those properties 
which depend on a restricted portion of the figure. For 
example, the total curvature of a surface at a point only 
depends on the shape of the surface at that point. Suc¬ 
cinctly we may say that Differential Geometry is the 
study of geometry in the small. This chapter is not intend¬ 
ed to be a complete course on the subject. However, 
sufficient theory is developed to indicate the scope and 
power of the tensor method. 

§ 88 . Permutation tensors 

In the Euclidean three-dimensional space let us in¬ 
troduce the quantities defined by 

( 88 . 1 ) e ilk — y/g e ijk ; e i,k = e ijk , 

where e tik are the permutation symbols defined by the 
following conditions :- 

(i) e tik = 0 if any two of the indices *,7 andfc are equal, 

(ii) e i23 = e 231 = e 312 = + 

(iii) e 18a = e 321 = e ai3 = 1 , 
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and g is the determinant formed from the fundamental 
tensor g (i of the space referred to some general coordinate 
system, which is not necessarily rectangular cartesian. 
The definitions show us that e iik , e m and e iik are skew- 
symmetric in all their indices. (Note that in this chapter 
the range of Latin indices is now from 1 to 8). 

We shall first prove that although e iik is not a tensor, 
both e m and e tik are tensors. We observe that 

dx* dx* dx k _ dx* dx* dx k _ dx* dx i dx k 

e * ,k d£ l dx m dS* e,i1c dx 1 d&" d£" giih dSP* dx 1 dx" * 


Thus e ijk 


dx*dx* dx k 


is skew-symmetric in l and m. Simi- 


dx l dx™ dx" 

larly it is skew-symmetric in the indices l, m and n. But 
this expression apart from sign is the Jacobian determi- 
Jljf J 

It therefore follows from the theory of deter- 


n “*i» 

minants that 


dx* dx* dx k __ 

6iik d£'dsrd&~ e ' mn 


dx T 


dx* 


Now the fundamental tensor g tf transforms to g i} when 
we transform to the coordinate system x*. We find from 
(21.1) that their determinants satisfy the equation 


i=S 


dx r 

dx* 


The quantities e iik transform to e ijk where 

__ y- dx* dx* dx k 
- Vg**m dfiftpidpi- 

Hence 


®Imn V? ®lnm V? 


Imn 


oaf 
dx* 


( 88 . 2 ) 


_ dx* dx* dx k 
g,m " - e<ik d&d&*d&’ 
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from which we see that e (ik is a covariant tensor of the 
third order. Also 

dx T _ 1 9® 1 dx m 9 ®". 

e,ron “ VI SF 0# Si 5; 

therefore, 


nlmn 


— « — 


VI 


i 

VI 


glmn 


_ ifk 9®* 9®" 9® n 

0®* 0®* 05* 


which shows that e lmn is a contravariant tensor of the 
third order. We call e m and e iik the permutation 
tensors. 

If the coordinate system is rectangular cartesian, g = 1 
and the permutation tensors have as their components 
the permutation symbols. In this coordinate system the 
covariant derivatives e {fk , and e m tl , which are both 
tensors, are zero. Thus the covariant derivatives e {ik> t 
and B l,k , are zero tensors in all coordinate systems. 
Hence the permutation tensors behave like constants 
with respect to covariant differentiation. 

From a covariant vector A t we can form the vector 

B* = e**'A ltk 

and we call B i the curl of the vector A { and write it 
curl A { . 


Ex. 1. Prove directly that = 0. 
Ex. 2. Prove that = gag /a gtne ,mn - 


Ex. 8. Show that the components of curl A ( are 
1 CbA t _ J_ IlA, _ _1_ [bAi _ ^£,\ 

Vi lx>) ’ Vi v>** to 1 ) * Vi V >* 1 *** /* 


§ 89. Vector product 
We can form the contravariant vector 
(89.1) C l = e^ k A J B k 

from the two covariant vectors A f and B { . To find the 
geometrical interpretation of the vector C*, let us choose 
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a rectangular cartesian coordinate system. In such a 
system the distinction between contravariance and 
covariance disappears and the C* are now the components 
of the vector product* of A t and B t . Thus C { is a vector of 
magnitude AB sin 6, orthogonal to both vectors A { and 
B it where 0 is the angle between these vectors. Its direc¬ 
tion is uniquely determined by the fact that A it B { and 
C* form a right-handed system. 


§ 40. Frenet formulae 

In this section we shall investigate the theory of 
twisted curves. Let the curve be given by the equations 
x i = £*(*), where the parameter s measures the arc- 
distance along the curve. Then the unit tangent vector 
T* to the curve is T* = dx*\ds and satisfies the equation 

gi/TiT* = 1. On differentiation we obtain g {i T J — 0. 

os 

That is, dTt/ds is a vector orthogonal to the tangent 
vector. Its magnitude k is called the curvature of the 
curve, and the unit vector 

(40.1) Ni = ~ir 

K 08 


is called the unit principal normal vector. We define 
the unit binormal vector B 4 to be the unit vector or¬ 
thogonal to both the tangent and principal normal vec¬ 
tors, and oriented so that the tangent, principal normal 
and binormal vectors form a right-handed system. Hence 
from (39.1) we have 

(40.2) B 4 = s iik TjN k . 

Since the vectors T { and N 4 are orthogonal, = 0. 

We differentiate intrinsically and obtain 


<3 N* dT 4 

«« r ‘-ar+S"i^' = 0 - 


* D. E. Rutherford, Vector Methods, p. 7. 
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dT* 

We substitute for from (40.1) and then replace 
gifNW 1 by both being equal to unity, and obtain 


e „T‘ (f? + «r<) = o. 


Now we differentiate g#NW # = 1, which gives us 
SN 1 

g ti N* — = 0. Thus 

g„ N i (|L' + ,3v) = 0. 

dN* 

Therefore the vector -=—(- kT* is orthogonal to both 

08 

the tangent vector T* and the principal normal vector 
N*. Accordingly it is in the direction of the binormal 
vector B*, and we can write 

(40.8) *' = 7 (^' + « T ‘) ■ 

The invariant t thus introduced* is called the torsion of 
the curve. Note that it may be positive or negative in 
contrast to the curvature which is essentially positive, 
being the magnitude of a vector. 

We differentiate (40.2) intrinsically and substitute 
from (40.1) and (40.8) and the result is 

— = 8 «fc ^N k + s«* T, —* 
ds 8 &8 SSk + B 8s 

= K T k ] 

which reduces on account of the skew-symmetry of e iik 
to 


8B* 

88 




But T it N { and B { form a right-handed system of unit 
* Some text books use r to denote the reciprocal of the torsion. 
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vectors; and so e iik T t B k =—N i and our equations become 

(40.4) ^ — tNK 

08 

Equations (40.1), (40.8) and (40.4) are called the 
Frenet formulae. On account of their importance in the 
theory of curves we group these formulae together for 
convenient reference in the form 


(40.5) 


6T*_ 
ds 
6N* 
ds 
dB { 
ds 


kN* 


— kT* + rB< 

— T N*. 


If the coordinate system is rectangular cartesian, the 
intrinsic derivatives become the ordinary derivatives and 
we have the well-known Frenet formulae. 


Ex. 1. Show that ic* = -. 

ds 8s 

AN 

Ex. 2. Prove that r =» s iik T { N S —- . 

' ds 

Ex. 8. A helix is defined to be a curve whose tangent vector makes 
a constant angle with a fixed direction. Prove that the 
necessary and sufficient condition that a curve be a helix 
is that the ratio of the torsion to the curvature is constant. 

Ex. 4. Prove that ‘ *Ii —* = . Hence a curve 

&8 Js 8 W ds \ K ) 

is a helix if and only if c«‘ — —i —* = 0. 

8s 8s* 8s* 


§ 41. Surface - First fundamental form 
We will now investigate the theory of surfaces. The 
three equations 

(41.1) x l = u 8 ); afi = x 2 (u\ tt 2 ); a 3 = <r*(tt l , u 8 ), 
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where u 1 and v? are parameters and the x i are three func¬ 
tions of u 1 and u 3 which are real and continuous, generally 
represent a surface. These equations may be briefly 
written x { = #*(«*) if we adopt the convention that Greek 
indices will always have the range 1, 2. A point at which 

[ dx*! . 

^ I is of rank two is called a regular 

point, A point may fail to be regular either because it 
is a singularity on the surface, for example, the vertex 
of a cone, or because it is a singularity of the parametric 
representation, for example, the poles of a sphere. (See 
equation (41.2)). 

Each pair of values of u a determines a point on the 
surface. That is, the u a form a coordinate system upon it. 
Thus an equation of the type f(u\ u*) = 0 must define a 
curve. Although a unique point on the surface corresponds 
to a fixed pair of values of u a $ the converse is not neces¬ 
sarily true. This is illustrated by the equations 

(41.2) x x = a sin^costt 2 ; aj^asinu 1 sinu 2 ; 0 s =a cos w 1 , 

which specify the surface of a sphere of radius a. In 
practice we restrict the range of the parameters so that 
the converse will also be valid. In our example, the 
parameters are limited to the ranges 0 u x ^ n and 
0 u 2 < 2tt. Then, except for the two poles, there 
corresponds to each point of the sphere a unique pair of 
values of u a . The equations u 1 = 0 and u l = n represent, 
not curves, but the poles respectively. At the poles, the 
coordinate u 2 is indeterminate. 

We shall exclude singular points from our discussion 
by considering only that portion of a surface on which 
there is a unique correspondence between its points and 
pairs of values of the coordinates u a . That is, in this 
portion every curve of the family u 1 = constant intersects 
every curve of the family u 3 = constant in one point 
only. The curves u 1 = constant are called u 3 -curves 
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and the curves u 2 = constant the u'-curves. Collec¬ 
tively they are designated the coordinate or parametric 
curves on the surface. Along a coordinate curve we 
choose the positive direction to correspond with increasing 
values of the variables u 1 or u 3 respectively. 

The contravariant vectors dx* and du* which represent 
the same displacement in space and on the surface respec¬ 
tively are connected by the equations 

(41.8) di r< = du* 

du* 


where we extend the summation convention to apply to 
Greek indices. Hence the line-element ds on the surface 
is given by 

ds 2 = g if das* dx* = g" ~ ~ du* dufi. 

Let us write 


(41.4) 


dx l dx* 

- g « dwfafi' 


from which it is clear that a a p is symmetric, and so we 
have 


(41.6) ds a = a a fdu*dufi. 

To this equation we apply the quotient law. Since a a p 
is symmetric it follows that a a p is a covariant tensor with 
respect to transformations of the coordinate system u*. 
We call it the fundamental surface tensor. Also 
a a pdu*dvfi is named the first fundamental form of the 
surface. 

Let us select rectangular cartesian coordinates in space. 
Then we readily obtain the metric in the well-known form 

ds 3 = Edu* + 2Fdudv + Gdv\ 
where u = u\ v = u 8 and 

E = Z — — F — £ — — r _ ydx* dx { 
i du du * i du dv ’ t dv dv 
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§ 42. Surface vectors 

When we transform coordinates in space, dx l is a 
contravariant vector but dw* is an invariant. On the other 
hand, if we transform the coordinates u a on the surface, 
dx i is an invariant but du a is a contravariant vector. 
Thus the equations (41.8) indicate that we can regard 
da^/du® as both a contravariant space vector and a co¬ 
variant surface vector. We may therefore introduce the 
notation 


(42.1) 


x' 


dx i 
d u a 


and rewrite (41.4) in the form 


(42.2) a afi = 


A curve on the surface is represented parametrically 
by the equations u a = u a (t). The vector du a jdt is a tangent 
vector to this curve. Its space components are then given 
by the equations 


(42.8) 


dx i dx { du a _ t du a 
dt d«“ dt ~ X * di ' 


But if the components of the space tangent vector dx i fdt 
are fixed, (42.8) consists of three equations for the two 
unknowns du a /dt. They are not consistent unless the vec¬ 
tor lies in the surface, when a unique solution would exist. 

Next we consider a surface vector-field A a . We can 
set up a unique curve C on the surface by the differential 
equations du a /dl = A a , provided that the vector-field is 
fixed at some particular point. Then A a is the surface 
tangent vector to this curve C. Let us designate the space 
components of A a by A 1 , and so equations (42.8) state 
that these components are connected by the relations 

A* = w^A a . 


( 42 . 4 ) 
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The magnitude of the vector A* is given by 
{AY = gii A*Ai = g if aia$A*Bft. 

That is, 

(42.5) (^)a = a afi A*At. 

In particular, du a Jds is the unit tangent vector to a curve 
on the surface if the parameter s measures arc-distance 
along it. 

The angle 6 between the two unit vectors A* and B { 
is obtained from 

cos 0 = g ii A i B ) = g {i co i a ai , fi A*B0. 

That is, 

(42.6) cos 6 = a a fiA a B0. 

It follows that the necessary and sufficient condition for 
the orthogonality of two surface vectors A a and BP is 

(42.7) a afi A*BP = 0. 

We see from equations (42.5), (42.6) and (42.7) that 
the familiar formulae apply equally well on the surface 
provided that we employ in them the components of the 
vector in the surface and the surface fundamental tensor. 
Also we can raise and lower indices of surface tensors in 
the usual way with the surface fundamental tensor a a p 
and its conjugate symmetric tensor a a P. These two tensors 
are connected by the equations 

(42.8) a aft a*y = d}, 

where <SjJ is the two-dimensional Kronecker delta. It is 
worthy of note that 

(42.9) a 11 = a M /a; a 12 = a 21 = — a 28 = a^Ja, 

where 

a = °ll a 22 (du) 8 - 
Ex. Prove that A a = x a T A r . 



68 TENSOR CALCULUS $ 43 

§ 48. Permutation surface tensor 
In section 88 we introduced the permutation tensors 
in space. Similarly we introduce on the surface the 
quantities defined by 

(48.1) e a 0 — yfa e afi ; e*0 = -^=. e afi , 

where 

% = ^22 = ^12 = " I " 1 } = - 1 . 

It is left as an exercise for the reader to show that e a p 
and are a covariant and a contravariant skew-symme¬ 
tric surface tensor respectively. They are called the 
surface permutation tensors. We see that e a fi can 
be obtained from by lowering the indices since 

Safi = a ay &fiS 

Now we derive an important formula for the angle 6 
between two unit vectors A a and B a . This angle is given 
by cos 6 = a a pA a BP. Accordingly, using (42.6) 

sin 2 0 = 1 — aafiA*Bfiay t Ay B s 

= { a ayO,fi» — a a fiOyt)A a Ay BP B 8 
— ae^i e 7 pA a Ay B? B s 
= e aa e y fiA*ArB0B s 
= (SosA'B*)*. 


To remove the ambiguity in sign the convention is made 
that we choose that value of 0 which satisfies 

(48.2) sin 0 = + e a tA a B t . 

In accordance with the convention just made, we say 
that the rotation from C a to D a is positive if the invariant 
e a pC a D 0 is positive. This in effect chooses the positive 
rotation as that one which rotates C a to D* through an 
angle less than or equal to n. 
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Let us form the contravariant vector 
(48.8) B» = 

from the covariant unit vector A a . Its magnitude is 
given by 

(B)“ — a a pB a BP = a a pev a e? fi A y At 

= + 2fl 12 e« e M^ i ^ 2 + o 11 (e“)»(^ 8 )* 

= {°2aMi) a 2a is^i^2 "I" ^liMa) 8 }* 

That is, in virtue of (42.9) 

( B ) s = o w (^i) 8 + 2a u A 1 A i + a™(A z )* 

= a^A a A fi = {A)* = 1. 

Thus B a is a unit vector. Further the angle 6 between 
A a and B a satisfies 


sin 6 = e afi A*BI> — e a ^A a B 0 = BPB fi = 1. 

Therefore 6 = jr/2 and by the above convention, equation 
(48.8) determines the unit vector B a orthogonal to the 
unit vector A a and oriented so that the rotation from 
A a to B a is positive. 

We now apply (48.2) to calculate the angle a> between 
the coordinate curves. The unit tangent vectors to the 

u 1 - and « a -curves are * <5? and —jL= <3? respectively. 

Hence the angle co satisfies 


(48.4) sin co = 


V<hi a , 




22 


\Zojja, 


°12 


22 


= t. 


a 


a n a 2* 


We also note from this equation that the rotation from 
the direction of a u 1 - curve to a u 2 -curve is always positive. 
It is easy to deduce that the necessary and sufficient 
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condition that the coordinate curves are orthogonal 
everywhere on the surface is that Oj a vanishes everywhere. 
We then say that the coordinates are orthogonal 
curvilinear. 


§ 44. Surface covariant differentiation 
We can form Christoffel symbols starting with the 
fundamental surface tensor a a p, and hence introduce 
covariant and intrinsic derivatives which are now sur¬ 
face tensors. So A a> p will denote the covariant derivative 
of A a with respect to uP, Written in full 



There will be no confusion between the Christoffel sym¬ 
bols formed with the g {i and those formed with the tensor 
a a fi since the Latin and Greek indices will clearly dis¬ 
tinguish which symbol is meant. We can also form the 
Riemann-Christoffel surface tensor from a a f and the 

curvature surface tensor R e fiyj = R% y s • 

Following the argument of section 26, we can show 
that the geodesics on the surface are the solutions of the 
differential equations 


(44.1) 


S (du a \ _ dPu a I a 1 duP du7 

Ts = ~W + fez dT ~dT 


0 . 


Corresponding to (27.1) the geodesics also satisfy the 
equation 


(44.2) 


du a duP 
a * fi ~ds ~ds “ 


Thus, in practice, we need only consider one of the two 
equations (44.1) together with the first order equation 
(44.2). 

As in section 28 we can introduce a system of geodesic 
coordinates on the surface so that at any particular point, 
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called the pole, all the surface Christoffel symbols are 
zero. At the pole, covariant and intrinsic derivatives 
reduce respectively to the corresponding partial and total 
derivatives. 

Again the theory of parallelism, outlined in section 29, 
applies to surface vectors. The vector-field A a is said 
to be parallel along the curve u a = u a (t) if 



Note that the covariant derivatives of a a fi, a*t and <5^ 
are zero. We cannot apply the method of section 88 to 
prove that E a p >r and are both zero because it is general¬ 
ly impossible to choose a rectangular cartesian system of 
coordinates on an arbitrary surface. Instead we select a 
geodesic coordinate system. At its pole, the Christoffel 
symbols are zero. Hence the partial derivatives of 
are also zero. Thus the partial derivatives of the deter¬ 
minant a are zero at the pole. That is e a ^, r and are 
both zero there. It follows immediately that these tensors 
are zero at every point on the surface and in all coor¬ 
dinate systems. Hence the permutation surface tensors 
behave like constants with respect to surface covariant 
and intrinsic differentiation. 

Ex. Show that the conditions that the u'-curvea and the u*-curves 

be geodesics are j^J = 0 and = 0 respectively. Simplify 

these conditions if the coordinate system is orthogonal cur¬ 
vilinear. 

§ 45. Geodesic curvature 
We discussed in section 40 the theory of the twisted 
curve in space. To this we now add the corresponding 
theory applicable to twisted curves lying on surfaces. 
Let the curve be given by the equations u a = u a (s), 
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where a measures arc-distance along the curve. The unit 
surface tangent vector is 


(45.1) 



Since is a unit vector, we have a a fit a tP = 1 and intrinsic 

differentiation yields us = 0. This shows that dt^fds 

is a vector on the surface orthogonal to the tangent vector 
f*. Let us denote the unit vector in the direction of 6t a }6s 
by n a . Then 

(45.2) ^ = cm* 


where a is an invariant called the geodesic curvature 
of the curve. We also call n a the unit normal vector 
in the surface to the curve, and choose its direction so 
that the rotation from f“ to »“ is positive. Thus 

e a0 t^nP = 1 , 

and equations (45.2) determine a uniquely in sign as 
well as magnitude. It is now clear from (48.8) and sec¬ 
tion 17 that 

nP = + e a ^t a 

and 

tP — — e a 0n a . 

Intrinsic differentiation of the first equation gives us 

~I~ = ~T = ae ?^ n * — — 

08 68 


Combining this with (45.2) we have the surface Frenet 
equations of a curve 


(45.8) 


6t? 

6s 


= an a 


6n* 

~8a 


= — ot*. 


Although i a is the same vector as T l , their components 
being connected by the relation T i = x i a V l , it is important 
to note that n a is in general neither the principal normal 
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N* nor the binormal B { . It does lie however in the normal 
plane of the curve determined by N { and B*. 

Along a geodesic of the surface di a /ds = 0 and hence 
o = 0. Conversely if o = 0, then dtP/ds = 0 and so the 
curve is a geodesic. Therefore, the necessary and sufficient 
condition that a curve be a geodesic is that the geodesic 
curvature be zero. 

Ex. Prove that the geodesic curvatures of the coordinate curves are 


§ 46. Normal vector 


We shall now derive an expression for the unit normal 
vector at any point on a surface. Let us choose its 
orientation so that the M 1 -curve, the u 2 -curve and the 
normal at the point form a right-handed system. The 
surface unit tangent vectors to the coordinate curves are 


1 1 

— <5J and —==. b\ respectively, and the corresponding 
V®11 V fl 22 l 1 

space components are _ a* and _ x l a 6 f, that is 

X i V°ii V^2a 

;— aj and ■ -p=. xi. In virtue of (48.4) and the co- 
V^u 2 V 

variant form of (89.1) we have 


That is, 
(46.1) 







It is not dear that £ { is a covariant vector because 
appears in this equation. However the vector form is 
clearly seen from the equivalent equation 

(46.2) = 
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For purposes of calculation (46.1) is the more suitable 
and can be readily written as a determinant, but (46.2) 
is preferable for theoretical purposes. Since the vector 
does not lie in the surface, there is no corresponding surface 
vector | a . We shall often require the important equation 

(46.8) g„S*xj, = 0, 

which expresses the fact that the normal is orthogonal 
to the surface vector This equation is also an im¬ 
mediate deduction from (46.1). 

Ex. The tangent surface of a curve is defined to be the surface 
generated by all its tangent lines, called the rulings. Show that 
the normals to the surface along a ruling are all parallel to the 
binormal of the curve at the point of contact of the ruling. 

§ 47. Tensor derivatives of tensors 
In the theory of surfaces we require tensors which 
possess both Latin and Greek indices, for example x* a . All 
such tensors in this chapter will be contravariant with 
respect to space but covariant with respect to the surface. 
We can then select A* a as a typical tensor. When we 
change the coordinate systems both in space and on the 
surface the transformation law is 

A< = A'— — 

* * dx* du a ’ 

We now ask what tensors can he constructed by differen¬ 
tiation? Here we follow section 80 closely. As our space is 
Euclidean the concept of parallelism in it does not depend 
on the selection of a curve. Accordingly choose an ar¬ 
bitrary parallel vector-field X { in space. On the surface 
take an arbitrary parallel vector-field Y a along the curve 
C whose parameter is t. Then 
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and 


tiY* 

dt 


+ 


f;,K£ r =°- 


Form the invariant A#X ( Y a . By differentiation and the 
use of the equations of parallelism we have 

+ M*) Xl % r °- Aix ‘{ey) Y '%- 

That is, on making an appropriate change of dummy 
indices, 

z lA - x ‘ Y ’ ] - K?+U) 4A ~ U) 4] t x - F *- 


Applying the quotient law we see that the expression in 
square brackets is a tensor, which we call the tensor 
derivative of A* a with respect to and we denote it 
by the semi-colon notation 

It is possible to choose a rectangular cartesian system 
in space and a geodesic system on the surface. Then at 
the pole, the tensor derivatives are the partial derivatives. 
Consequently the laws of tensor differentiation are the 
same as those that apply to covariant derivatives. 

We must now extend the concept of tensor differ¬ 
entiation further to space tensors and to surface 
tensors. It is clear that the tensor derivatives of surface 
tensors are identical with their covariant derivatives. 
Following the method of this section, we see that the 
tensor derivative of a space tensor with respect to u® is 
the tensor obtained by inner multiplication of its 
covariant derivative with respect to x l by the tensor a£. 
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For example, A i! . a = A {i k x%. Thus the tensor deriva¬ 
tives of g,„ g«, 6 { f , e ijk , e iik , a afi , a*t, d%, e afi and are 
all zero. That is, they may be treated as constants with 
respect to tensor differentiation. 


§ 48. Second fundamental form 
By tensor differentiation 
d z x i 


(48.1) 


duf*du0 


+ 


f * I X 1 xl — f y \ X* 

\jk] X '* x 0 w\ y> 


which shows that x { a .j is symmetric in a and That is, 
:^. a . Now the tensor derivation of (42.2) yields 




= x) 


= °- 

We subtract this equation from the sum of the two 
similar equations obtained by cyclic interchange of a, /? 
and y, Because is symmetric the result is 

(48.2) = °- 

This shows that ** ;>J is a contra variant space vector 
orthogonal to all vectors lying on the surface. Thus 
it is co-directional with the normal vector Hence 
quantities b a p must exist so that 

(48.8) 

Further, it follows that b a ^ form the components of a 
covariant symmetric surface tensor. Equations (48.8) 
are known as Gauss’s formulae. Since £* is a unit 
vector, inner multiplication of (48.8) by and sub¬ 
stitution from (46.2) yields 


(48.4) b a p 2' e ^ e iik X a;f) X y X » — y/~Q, • 

The quadratic form 

(48.5) bafiduPdvtA 
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is called the second fundamental form of the sur¬ 
face. It is now possible to construct the invariant 
(48.6) H = 

which is called the mean curvature of the surface. 

Ex. If the space coordinates are rectangular Car tesian , show 
that 

. 1 h ** 1 t * 

§ 49. Third fundamental form 
Tensor differentiation of the identity guPP — l yields 
“S — °. That is, is a contravariant space 

vector orthogonal to the normal vector. Accordingly it 
lies in the surface. Hence quantities r)£ exist so that 

( 49 - 1 ) £*ia = yilxy 

The quotient law then states that rj£ form a mixed surface 
tensor. We now differentiate (46.8) tensorially and obtain 

+ giiPVfiia — °> 

which reduces on substitution from (49.1) and (48.8) to 

giiVWy^ + giiPbafiP = 0 . 

We apply (42.2) and the result is 
(49.2) b afi = — a fiy if a . 

Inner multiplication by aP a yields 

( 49 *®) Va~ — aP*b a fl ; 

and so we can rewrite (49.1) in the form 
(49.4) = — ayfi b ar x { fi . 

These equations are known as Weingarten’s formulae. 
Introduce the symmetric surface tensor 

( 49 - 5 ) c *fi = g<i£a£U 
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and form from it the quadratic form c a pduf*dufi which is 
then called the third fundamental form of the surface. 

Ex. Prove that c a fi — a?* b a y bpg . 

§ 50. Gauss-Codazzl equations 
We are now in a position to obtain the central formulae 
in the theory of surfaces. For the moment, choose the 
space coordinates to be rectangular cartesian and the 
surface coordinates to be geodesic. Then at the pole, 
tensor differentiation of (48.1) yields 

t _ d 3 x* _ d_ | a) f 

~ dufidufidu? du? \a/?J 

Thus 



The expression in square brackets is the surface Riemann- 
Christoffel tensor at the pole. Therefore we have set 
up a tensor equation 

(50.1) <* y — < yj » = R'.fX 

which must be true at every point of the surface and in 
all coordinate systems. 

We substitute from (49.4) in the tensor derivative of 

(48.8) and obtain 

x a;0y = ^aff; 0, ea b a fib ye x r' 

We can consequently write (50.1) in the form 
(Pafl; y — — ° e<r (^a£^y« &^e)®Jr = R?afly x ir’ 

Inner multiplication by £< and in virtue of (42.2) 

and (46.8) yield the respective equations 

(50.2) ^afi;y &ay;ft ~ 0 
and 

(50.8) 


Rpafiy — bay bpff b a fl bpy- 
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The reader is asked to verify, that (50.2) consists of 
only two independent partial differential equations. They 
are called the Godazzi equations. As there is only one 
distinct component of the curvature tensor in two dimen¬ 
sions, (50.8) reduces to the single equation 

(50.4) *ms = M2a— (&ia) 8 

which is called the Gauss equation. By means of (85.2) 
we can write this equation 

(50.5) K = — 

a 

where b is the determinant formed from the b a p and K 
is the Riemannian curvature of the surface. On a surface 
it is more usual to call K the Gaussian or total cur¬ 
vature. 

It can be proved* that a surface is uniquely deter¬ 
mined except for a translation or rotation in space when 
the first and second fundamental forms are given. This 
theorem can be precisely formulated as follows;- If a a p 
and bap are given functions of u 1 and « a , then there exists a 
surface x* = x*(u a ), uniquely determined except for its 
position in space, which has a a ftdu a dvfi and b a pdu a duP as 
its first and second fundamental forms respectively, 
provided that a a pdu a duP is positive-definite and that a a p 
and b a ft satisfy the Gauss-Codazzi equations. 

Ex. Prove that e a p = 2 Hb a fi — Ka a 0 • 


§ 51. Normal curvature-asymptotic lines 
On a surface consider the curve u a = «“(«), where s 
measures arc-distance. The equations of the curve in 
space will be x i = x i (s). Then the space vectors T i , N* 
and B* and the surface vectors tf* and «“ at any point of 
the curve satisfy the Frenet formulae (40.5) and (45.8). 

* L. P. Eisenhart, Differential Geometry, p. 157—159. 
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The tangent vectors T l and.#* are connected by T l = 
Intrinsic differentiation yields in virtue of (48.8) 


— = T l — = T i x>— = T i — 
ds •* ds •* * ds ds 

* ds aifi ds 


Applying the Frenet formulae we have 
kN* — an a Xa + b a 0t a tft 


which becomes 


(51.1) kN* = an 1 + b afi t* tfi? 


when we designate the space components of n a by n*. 
Let us introduce the angle 0 between the principal normal 
N* and the surface normal Then inner multiplication 
of (51.1) by £, yields 

(51.2) /ccosfl = baft*#, 


since $ { is orthogonal to the vector n* which lies in the 
surface. The invariant b a pi a tfi is the same for all curves 
which have the same tangent vector if* at the point on 
the surface. Accordingly we deduce Meusnier’s theorem 
‘For all curves on a surface which have the same tangent 
vector, the quantity k cos 6 is constant’. This quantity 
is called the normal curvature at the point and is 
denoted by «(„). Hence 


(51.8) 


K(n) = bafiPtfi = 


b a pdu*dufi 
a a pdu a dvfi ’ 


If we choose a plane section through the normal to the 
surface, then 0 = 0 or n. That is «•(„) = ± k. So the 
normal curvature in any direction is equal in magnitude 
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to the curvature of the normal plane section of the surface 
in that direction. 

Along a geodesic a = 0, therefore equations (51.1) 
reduce to kN* = b afi ^tfi^. Thus either k = 0 or N { = 

We deduce that a geodesic on a surface is either a straight 
line or is a curve whose principal normal is co-directional 
with the surface normal at every point. Conversely, if 
N* = ± then inner multiplication of (51.1) by n ( 
yields a = 0. That is, the curve is a geodesic. 

The directions at a point on a surface which satisfy the 
equation b a fidu a du0 = 0 are called the asymptotic 
directions. If all the points of a curve have asymptotic 
tangent directions, the curve is called an asymptotic 
line. The asymptotic lines on a surface are given by 
b a pdu a du0 = 0. Thus along an asymptotic line we have 
kN* = an*. Hence, since N* and n* are both unit vectors, 
either k = a = 0 or the curvature and geodesic curvature 
of an asymptotic line are equal in magnitude and its 
principal normal lies in the surface. Consequently the 
binormal of an asymptotic line, which is not a straight 
line, is co-directional with the surface normal. The con¬ 
verse is also valid. 

Ex. Prove Enneper’s formula tlxat the torsion t of an asymptotic 

line is ± V — K where £ is the Gaussian curvature of the 

surface. 

§ 52. Principal curvatures - lines of curvature 

The normal curvature /c (n) of a surface in the direction 
t* is given by /C(„) = b a pt a t0, where the tangent vector 
satisfies a^t a tP = 1. The maximum and minimum values 
of «(„) can then be determined by the method used in 
section 19. They correspond to the principal directions 
determined by fe a/ j and are given by the roots of the deter- 
minantal equation | b a p — Xa a p | = 0, which reduces in 
virtue of (50.5) and (48.6) to 

(52.1) A a — 2HX+K = 0. 
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The roots /C(u and k$) of this equation are called the 
principal curvatures of the surface at the point. The 
principal directions and Ufa corresponding to the prin¬ 
cipal curvatures at the point respectively satisfy 

(b a fl i) “ 0, 

{bafl *( 2 ) a afi)tfz) 5=5 

A point at which K(d = ac( 2 ) is called an umbilic. At 
all other points section 19 tells us that and are 
orthogonal to one another. A curve on the surface which 
is a principal direction at all of its points is called a 
line of curvature. 

At an umbilic, equation (52.1) has coincident roots. 
That is fl 2 = K and the reader may verify* that this 
result can be written 

4a(o 11 

+ [^11(^1 ^22 — a 22 ^ n ) — ^s^u)] 2 ® 

Since a a pduPdufi is positive-definite, a is positive, and we 
deduce that 

^ 1^12 — ^12 ^11 = ** 11^22 a 22^11 ^ 

and so 

_ ^12 _ ^22 * 

**22 

Equation (51.8) now shows that K( n ) is independent of 
the direction du a jds. That is, at an umbilic, the normal 
curvature is the same in every direction**. 

Ex. 1. Prove that the lines of curvature on a surface are given by 
eT* o ay du a dufi = 0. 

Ex. 2, If the coordinate curves are lines of curvature, prove that 
a xl = b lt « 0, and conversely. 

Ex. 8, Prove that a surface, all of whose points are umbilics, is a 
sphere or a plane. 

* L. P. Eisenhart, Differential Geometry, p. 119. 

** The reader may find it interesting to relate some of the 
results of this chapter to chapter II of Rutherford’s Vector Methods. 
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§ 58. Orthogonal transformations 
In this chapter our aim is to present the Theory of 
Elasticity and so we restrict ourselves to a Eu clidean 
space of three dimensions. In it we choose a right-handed 
system of rectangular cartesian coordinates and we denote 
them by y v where Latin indices have the range 1 to 8. 
The line-element ds is then given by 

(58.1) ds 2 = dy { dy { = 8 {j dy t dy f , 

where 8 if is the Kronecker delta. 

The linear equations, (compare (6.8)), 

( 58 - 2 ) yi = auyi + b t , 

where b { form a set of three constants and a {) a set of 
nine constants define a transformation to a new coor¬ 
dinate system . The necessary and sufficient conditions 
that the y t form a set of rectangular cartesian coor¬ 
dinates is 

d * 8 = dy t dy { = a u a {k dy,dy k = 8 jk dy t dy k . 

That is, 

(«««« — 8jic)dy,dy k = 0 

for all values of dy t . Hence 
(58.8) = 8 f k . 

Inner multiplication of (58.2) by a ik yields the solution 
(58.4) y k = a {k §t — a ik bt . 
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Thus 

(58.5) 


%ji = dyi 

tyi &9i 


Qu¬ 


in virtue of these equations, we see by exa m in in g (9.1) 
that the distinction between contravariance and co- 
variance has disappeared. Accordingly we shall write all 
indices as subscripts on condition that we allow only 
transformations of the type (58.2) subject to (58.8). We 
have already anticipated this in our notation y ti 6 {j and 
o w . However we may sometimes wish to adopt a cur¬ 
vilinear coordinate system such as spherical polars. It is 
then necessary to reintroduce the distinction between 
contravariance and covariance. This will be indicated by a 
return to the coordinates x*. 

The transformation (58.2) is equivalent to the com¬ 
bination of the two transformations y { = y\ + b t and 
y\ — a aVi- The transformation i} f — y' t -f- b { merely 
defines a translation to new parallel axes. The transfor¬ 
mation 


(58.6) y\ = a (i y, 

subject to the six conditions (58.8) is said to define an 
orthogonal transformation. It follows from (58.8) that 
the determinant | Oy | of an orthogonal transformation 
is either +1 or — 1 and we say that (58.6) defines a 
positive or a negative orthogonal transformation respec¬ 
tively. It is well known* that the y\ axes form a right- 
handed or a left-handed system if the orthogonal trans¬ 
formation is positive or negative respectively. Further a 
positive orthogonal transformation defines a rotation of 
the axes about the origin. 

An alternative set of equations to (58.8) is obtained 
by considering 

* W. H. McCrea, 'Analytical Geometry of Three Dimensions’, 

pp. 10—18. 
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ds* = dy i dy i = = d^dy^y^ 

from which we deduce that 
(58.7) a H a lci = &ik’ 

Ex. Show that in the rotation of axes defined by the positive 
orthogonal transformation (53.6), a u is the cosine of the angle 
between the y' { axis and the y s axis. The components a i{ are 
therefore the direction cosines of the y' { system with respect 
to the y f system. 

§ 54. Rotations 

Equations (58.2) may be interpreted from another point 
of view. We could say that they transform the point 
P whose coordinates are y { into the point P whose coor¬ 
dinates are y { referred to the same system of rectangular 
cartesian axes. We then call (53.2) an affine transfor¬ 
mation*. If in addition conditions (53.3) together with 
I a u I = 1 are imposed, we have the most general rigid 
body motion consisting of a rotation followed by a trans¬ 
lation. 

Next we wish to obtain the orthogonal transformation 
which represents a right-handed rotation through an 
angle about a line through the origin whose direction 

* Affine Geometry. According to Klein's Erlangen Programm, 
(F. Klein, Math. Ann. Vol 48, (1898), p. 68), a geometry comprises 
a system of definitions and theorems that express properties which 
are invariant with respect to a given group of transformations. For 
example, if the group of transformations consists of all rigid body 
motions, namely translations and rotations, then the geometry is 
termed metrical. This is the geometry of Euclid, less the similarity 
theorems, and its most important concepts are distance and angle. 

Let us now examine affine geometry, defined by the group of 
transformations (53.2), Suppose A i to be the vector joining the point 
with coordinates to the point with coordinates y t . Then 
A< = y { — z i9 and we immediately deduce from (58.2) that the 
transformation law of vectors is A t = a u A t . We now choose two 
parallel vectors A t and B i . The requisite conditions are 
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cosines are l { . By a right-handed rotation, we mean that 
a right-handed corkscrew will move outwards from the 
origin along the direction l { when it is rotated through 
an angle ip about this direction. Using the vector notation 
we have (fig.) 



OP - OP + PL + LP, 

where Q is the perpendicular from P on the line through 
the origin whose direction cosines are l i7 the angle 


AJB l = AJB 3 = AJB 3 . But each fraction equals ( a ll A l + a lt 
A t +a ll A z )l(a ll B l +a 13 B 3 +a li B z ) = AjB lf and similarly equals 
AJB Z and A 3 /B a . Thus we have AjB l = A^jB^ *= A 3 lB 3 . That 
is, the vectors A t and B { are parallel. Hence the parallelism of 
vectors is invariant in affine geometry. 

Further we see from (53.2) that every finite point transforms 
into a finite point. Therefore the plane at infinity is invariant, and 
so we can distinguish between a non-central and a central quadric 
according as the plane at infinity does or does not touch the quadric. 
However we cannot define either angle or distance in affine geo¬ 
metry, because they are not invariants under the affine group of 
transformations (58.2). 
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PQP —y, QP = QP, the plane PQP is orthogonal to 
l { and PL is the perpendicular from P to PQ. Since 

OP = y it we have OQ = lil k y k and thus PQ = l { l k y k —y { . 

Therefore ~PL - (1 —cos ip)PQ = (1 — cos y>)(l { l k y k —yi). 

Further LP is orthogonal to both PQ and l { in such a 

way that PQ, l { and LP form a right-handed system. 

Then by (89.1) the unit vector in the direction LP is 
given by e tik (l f l m y m — yj )l k /PQ = — e ( , kVi l k IPQ. Hence 

LP = PQ sin ip {— e ijk y,l k /PQ) = e ijk l,y k sin ip. Thus we 
have 


Vi — Vi + (1 — cos y>){l { l k y k — y { ) 4- sin ip e iik l f y k , 
which can be written 
(54.i) y t = a ik y k 

where 


(54.2) a ik = cos ip 8 ik -j- (1 — cos ip)l ( l k + sin ipe i)k l f . 


In the theory of elasticity we shall be particularly 
interested in infinitesimal rotations, in which case 
cos y — 1 and sin ip ^ ip. The infinitesimal rotation is then 
represented by 


That is, 

(54.8) 

where 


Vi = Vi + ve {jk l f y k . 
9t = y t + s ik y k , 


(54.4) 8{k = ipe ifk l f . 

It is clear that 8 tk is anti-symmetric. Conversely, let us 
consider equations (54,8) given that 8 ik is anti-symmetric. 
Then equations (54.4) comprise only three equations for 
the three unknowns l t , whose solution is in fact l t = — s^jip, 
h = — hi/V l z = — SxJip. Therefore equations (54.8) 
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always represent an infinitesimal rotation if ip is in¬ 
finitesimal and s ik is anti-symmetric. 

Ex. Calculate a a corresponding to a rotation of S0° about the y t 
axis. 


§ 55. Cartesian tensors 

A Cartesian tensor of the Af-th order in a three- 
dimensional Euclidean space is defined as a set of 8 W 
quantities which transform according to equations (9.1) 
when the coordinates undergo a positive orthogonal 
transformation. This is a less stringent condition than 
that imposed on a tensor. So we see that all tensors are 
Cartesian tensors but a Cartesian tensor is not necessarily 
a tensor in the usual sense. In virtue of (58.5) we have 
that A* t r. is a Cartesian tensor of the M -th order 
if the transformed components satisfy 

(55.1) •• k M’ 

on change of the coordinates by the positive orthogonal 
transformation y i = a ii y i . We see from (58.6) that both 
y i and their differentials dy f are Cartesian vectors. Also 
the Kronecker delta is a Cartesian tensor of the second 
order because 

— a ir a ii^T3 = a ir a ir ~ &ii 

in virtue of (58.7). Further we deduce from (88.2) that 
the permutation symbol e ijk is a Cartesian tensor of the 
third order. Thus a ik and s ik introduced in the last 
section are Cartesian tensors of the second order. 

The fundamental tensor of the Euclidean space is the 
Kronecker delta 5 ti . Hence all the Christoffel symbols are 
zero and so the comma notation for covariant derivatives 
now denotes the familiar partial derivatives which are 
Cartesian tensors. 

The use of (55.1) instead of (9.1) shows us clearly that 
the quotient law of section 12 still applies to Cartesian 
tensors. 
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§ 56. Infinitesimal strain 

Consider a body which is strained by the action of 
applied forces. The particle at the point P with coor¬ 
dinates y t referred to a rectangular cartesian system is 
displaced to the point P with coordinates y { + it,. Similar¬ 
ly the particle which was at Q with coordinates z { is 
displaced to the point Q with coordinates z { -f- v t . We 
define the extension &(pq) of the straight line joining the 
unstrained points P and Q to be the increase in length 
per unit of length due to the strain. That is 


(56.1) 
We have 

and 


e {PQ ) = 


PQ-PQ _PQ 
PQ PQ 


1. 


(PQ) 2 = (Vi — z { )(y { — z { ), 


(PQ ) 2 = (Vi + «< — — v { )(y { -f u { — Zf — v t ) 

— (Vi — *i)(y { — z { ) + 2(y { — zJiUi — v t ) 

+ («< — »<)(«< — 
= (P@) 8 Jl — 2l i( U i~ V i) _L («* — *>«)(«« — *<)\ 


PQ 


(PQ) 2 


where l ( = (z^ — Vi)\PQ are the direction cosines of the 
unstrained line PQ, Hence 


e (PQ) = 


2lj(u t — V i) t (u t — »,)(«,— 

PQ (PQ) 2 / 


We confine our attention to the case of infinitesimal 
strain, where the assumption is made that the displace¬ 
ment vectors u, t and v { are small compared to PQ. 
Neglecting quantities of higher order than the first, we 
obtain 


e (PQ)= ~ 


li(u< ~ Vi) 

PQ 


Let us now choose Q to be in the neighbourhood of 
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P so that y t — z { is small. Then by Taylor’s theorem for 
a function of three variables* 


v { — u t + ( Zf — + terms of higher order in z t — y } . 

Neglecting terms of higher order than the first in z f — y if 
we obtain that the extension e at P in the direction 
determined by the unit vector l { is given by 


e = 


~ Vi) 

PQ 


u t,i — 


Introduce the symmetric Cartesian strain tensor 
e if by the equations 


(56.2) e {l = + t* M ), 

(there is no confusion between e { f and the permutation 
symbols of the previous chapter which are distin¬ 
guished by Greek indices), and we finally obtain the 
extension e as the quadratic form 


(56.8) e = 

The dilatation or expansion 0 is defined as the 
increase in volume per unit of volume. That is 


(56.4) 


AV — AV AV 
~ AV ~AV 


where AV denotes the strained volume corresponding to 
the volume AV. But 


AV 

AV 


%< + «<) 


d(y,) 


= 1 + “i.i + « a .a + «3.s + terms of 

higher order. 


Thus in virtue of (56.2) and (56.4) we have 
(56.5) 0 = e it . 


The components of the strain tensor are not entirely 
arbitrary. To prove this, differentiate (56.2) twice and 
obtain 


* R. P. Gillespie, Partial Differentiation, p. SO. 
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e u,ki = 

from which it immediately follows that 

(56.6) e if'kl + e kl,ii — e ik,H — e jl,ik — 0. 

There appear to be 81 of these compatibility equations. 
In actual fact some are repeated due to the symmetry of 
the strain tensor and others are satisfied identically. The 
reader is asked to verify that only six of these equations 
are independent. 

To conclude this section, we investigate several im¬ 
portant examples of strain. 

(1) Uniform dilatation. Consider the displacement 
vector u { = cy { where c is a constant. We have e (i = cd ti 
and the dilatation 6 = Sc. Thus the extension at any 
point in any direction is constant and equal to one-third 
of the dilatation. 

(2) Simple extension. Consider the displacement 
vector Ui — clilfyf where c is a constant and l { a unit 
vector. We have u i f = cl { lj and so e i} = cl t l t and the 
dilatation 6 = cl { l { = c. Also e ti l t l f = c and thus there is 
an extension at any point in the direction l t of amount 
equal to the dilatation. The extension in any direction 
orthogonal to l t is easily seen to be zero. If c is negative 
we refer to it as simple contraction. 

(8) Shearing strain. Consider the displacement 
vector u i = 2cl i m f y j where c is constant and both l { 
and are unit vectors. A brief calculation yields 
e u — c {h m j + h m t) an d the dilatation 6 = 2cl i m { . Conse¬ 
quently the dilatation is zero if the directions l { and m { 
are orthogonal. 

Ex. Show that a simple extension along any direction together with 
an equal simple contraction along an orthogonal direction is 
equivalent to a shearing strain along a direction bisecting the 
angle between the given directions. 
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§ 57. Stress 

The forces acting on a body are either external or 
internal. The external forces may consist either of body 
forces such as gravity which act on every particle of it, 
or of surface forces which act on the external surface of 
the body, for example the pressure between two bodies 
in contact. If F f denotes the body force vector per unit 
volume, then the force acting on an element of volume 
AV is FfAV. Similarly if T { denotes the surface force 
vector per unit area, then the force acting on an element 
of surface AS is T t AS. In order to discuss the internal 
forces, we select a small element of area AS inside the 
body and denote the direction cosines of the normal to 
this element, which is approximately planar, by »*. We 
call one side of the element AS positive and the other 
side negative. Then the action of the positive side on the 
negative side is the internal surface force T } AS where 
T f is the force per unit area on the element AS. It is 
called the stress vector and is in general a function of 
the coordinates of the point which determines the position 
of the element AS and of the direction cosines n f of the 
normal to AS. It is well to emphasise that T t is not 
necessarily co-directional with n*. At all points on the 
external surface of the body T t becomes the external 
surface force. 

Consider a small rectangular parallelepiped with vertex 
at the point P whose edges are parallel to the coordinate 
axes. We form three stress vectors T^, T (a )i an d Tp) f 
corresponding to the small elements of area through P 
which are parallel to the coordinate planes. The stress 
vector T( t )j will be called positive if it acts in the positive 
direction of the y ( axis provided that the external normal 
is co-directional with the positive y { axis. If, however, 
the external normal is co-directional with the negative 
y { axis, then the stress vector T Wi is positive if it acts 
in the direction of the negative y t axis. In other words, 
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a stress which tends to stretch will be regarded as positive, 
whilst a stress which tends to compress is regarded as 
negative. We define nine quantities E {j by the equations 

(57.1) E ti = T Wf 

and we shall show that E tj is a Cartesian tensor, called 
the stress tensor. 

Construct the small tetrahedron PA 1 A i A a such that 
the edges PA { are parallel to the y t axes. The forces 
acting on the tetrahedron are the body force F f AV, the 
surface forces E it AS t (no summation over i) on the 
faces opposite to A { and the surface force T { AS acting 
on the face A 1 A a A a , where AV is the volume of the 
tetrahedron, AS t is the area of the face opposite to A { 
and AS is the area of the face A 1 A 2 A 3 . Let the positive 
direction of T t be that of the normal drawn outwards 
from the tetrahedron and whose direction cosines are n } 
and let p be the perpendicular distance from P to the 
face A 1 A 2 A 3 . Then AV = frpAS and AS ( = n ( AS. The 
equations of equilibrium obtained by resolving forces 
parallel to the axes are now 

F t AV— E ti AS { + T t AS = 0 , 

where the E it occurs with a negative sign, because the 
external normals are in the directions of the negative 
axes. We now substitute for AV and AS ( , then divide by 
AS and proceed to the limit as the tetrahedron shrinks 
to zero, in which case p tends to zero. The result is 

(57.2) T, = .E„n,. 

It follows from the quotient law that E if is a Cartesian 
tensor, and we can by (57.2) calculate the stress vector 
at any point, corresponding to any direction at that point 
in terms of its direction cosines and the stress tensor. 
We now cite several important cases of stress 
(1) Normal Stress. The vector T f is co-directional 
with ± We see from (57.2) that E if = Cd (j where C 
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is a constant. Hydrostatic pressure is an example of 
normal stress for which C is negative. 

(2) Simple Tension. Consider the stress tensor 
E if = Cl { lj where C is a constant and l { is a unit vector. 
Then the stress vector in the direction l { is T i =Cl i l i l i =Cl i 
and it is thus co-directional with ± h- However the stress 
vector in the direction m t orthogonal to Z< is T t = Cl, !,»»<= 0. 
If C is negative the stress is called a Simple Com¬ 
pression. 

(8) Shearing Stress. This is specified by the stress 
tensor E ({ = + l/tn<) where C is a constant and 

l t and m ( are unit vectors. 

Ex. Show that a simple tension along any direction together 
with an equal simple compression along an orthogonal direc¬ 
tion is equivalent to a shearing stress along a direction bisecting 
the angle between the given directions. 


§ 58. Equations of equilibrium 
Let us consider a body of volume V in equilibrium, 
which is enclosed by the surface S. We resolve forces 
parallel to the axes and obtain 

f v F t dV+f s T t dS- 0. 

Substitution from (57.2) yields 

j v F i dV+j s E ii n i dS = 0, 

which on application of Gauss’s theorem * becomes 

f r Pidr+f r B„',dv~ 0 . 

That is, 

J v {Ff + E ii t )dV = 0 . 

This equation is an identity, being true for any volume V. 
Hence 

(58.1) Ff + Eh' i = 0. 

* D. E. Rutherford, Vector Methods, p. 74. 
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The moments of a force about the axes are the com¬ 
ponents of the vector product of the force vector and the 
position vector of any point on the line of action of the 
force. Thus in tensor notation, the moments of the force 
F { about the axes are represented by the vector e ijk yjF k . 
Now take moments about the axes for our body in 
equilibrium and the result is 

j v e m y,F k dV + l s e {ik y,T k dS = 0. 

We substitute from (57.2), apply Gauss’s theorem and 
obtain 

j v e mViF k dV + j y (e ifk yjE lk ) t dV = 0. 

That is, 

j v e weVi(F k + E, kl )dV + j v e iik <5 #l E lk dV = 0. 

The first integral is zero in virtue of (58.1) and we have 

jy e iikFi k dV = 0 . 

This equation is also an identity and consequently the 
integrand e ilk E lk vanishes, from which we deduce that 

(58.2) Eu = E fi . 

Thus the stress tensor is symmetric and the elastic 
equilibrium equations for a body are given by (58.1). 

§ 59. Generalised Hooke’s law 
In the elementary theory of elasticity, Hooke’s law 
states that the tension of an elastic string is proportional 
to the extension. In other words, stress is proportional 
to strain. The corresponding assumption in the general 
theory of elasticity is that the stress tensor is a linear 
homogeneous function of the strain tensor. That is 

(59.1) E iS = c ijkl e kl . 
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It follows from the quotient law that c (fJel is a Cartesian 
tensor of the fourth order, and it is called the elasticity 
tensor. Further from the symmetry of E it and e kl we 
find that c im is symmetric not only with respect to the 
indices i and j but also with respect to k and l. 

A body is said to be homogeneous if the elastic pro¬ 
perties of the body are independent of the point under 
consideration. This means that the components of the 
elasticity tensor are all constants for a homogeneous body. 
We call a body isotropic if the elastic properties at a 
point are the same in all directions at that point. This 
means that the elasticity tensor c ijkl transforms to c im 
itself under any rotation of axes. 

§ 60 . Isotropic tensors 

A Cartesian tensor which transforms into itself under a 
rotation of axes is called an isotropic tensor. We have 
already met two isotropic tensors, namely d i{ and e iik . 
We shall now search for the most general isotropic 
tensor c im of the fourth order. Its transformation law 

(55.1) becomes 

( 60 . 1 ) c {jkl = a ir a jt a kt a lu c rstu . 

Let us rotate the axes through 180° about the y z axis. 
We deduce from (54.2) that a tk =— d« + 2 l t l k . But 
for this transformation l x = l 2 = 0 and l 3 = 1, and so the 
only non-zero components of a ik are 

= 1» % = If % = *4* 1* 

By direct substitution in (60.1) we obtain c im = — c iikl , 
that is, c im = 0 in the following cases 

(1) any three of the indices equal to 1 and the other 
equal to 8. 

(2) any three of the indices equal to 2 and the other 
equal to 8. 

(8) any two of the indices equal to 1, another equal to 
2 and the other equal to 8. 
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(4) any two of the indices equal to 2, another equal 
to 1 and the other equal to 8. 

Similar results are obtained by the corresponding 
rotations through 180° about the y 1 and the y z axes. 
Hence the only components which survive are those 
where the four indices are equal or equal in pairs. 

Let us now rotate the axes through 90° about the y 3 
axis. For such a transformation we deduce from (54.2) 
that a ik — l t l k -j- e m l } and so the only non-zero com¬ 
ponents of a ik are 

Oj2 = 1; Oji = + 1; Ojg = + 1. 

By direct substitution in (60.1) we discover that 

Cnn = C2222» 

c 1122 = c 22111 Cji33 = C2233, C33H = C3322, 

c 1212 — c 2121> Cm3 = C2328 > C3131 = C3232 > 

c 1281 = Cana, c 1331 = c 2338» C3113 = c 3228* 

Similar results are obtained by the corresponding rotations 
through 90° about the and the y z axes. Hence we can 
collect all our results in the form 


(60.2) 


c mt — c uu » 

C UH = C Ukk = c lli) — c llkk » 
C UU = c ikik — C UU = c lklk * 
c iiii — c ikki = c liil — c lkkl » 


where i, j, k and l are unequal and the summation con¬ 
vention does not apply. All other components are zero. 
The most general solution of equations (60.2) is then 


(60.8) Cff j.j — Xdijd n + fiducdjt + vd {l d fk + K&uki , 


where A, fi, v and k are Cartesian invariants, and & iikl = 1 
if all four indices are equal and otherwise zero. 

Finally we carry out a small rotation, which is re- 
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presented by a ik = d ik + s iJa where a it is skew-symmetric 
and of the first order in small quantities. Substituting 
in (60.1) and retaining only terms of the first order 
we have 

8 ir C rm + 8 i» c tikl + 8 ki C iUl + 8 lu C iiku = 

Select* = 2,j = k = l = l. Then in virtue of s 13 — — , 

the non-vanishing terms of this equation yield 

C 1U1 = C 22U H" c 2121 + c 2112 * 

Substitute in this from (60.3) and the result is 

X+p+v+K=X+p+v 
which gives k = 0, so we can now write 
(60.4) c ifkl = Ad {i d kl + fi8 {k d it + vd a 8 ik . 

This tensor is obviously isotropic and so it is the most 
general isotropic tensor of the fourth order. 

Ex. 1. Prove that the most general isotropic tensors of the second 
and third orders are X8 { , and respectively where X is 
an invariant. 

Ex. 2. Prove that c Uk , defined by (60.8) satisfies the symmetry 
relations Cmi — Quo and nit ■ 

§ 61 . Homogeneous and isotropic body 
In this section we confine our attention to a body 
which is both homogeneous and isotropic. Then the 
elasticity tensor is isotropic with constant components. 
We therefore require an isotropic tensor c iikt which is 
symmetric in both i and j and in k and l , and whose 
components are all constants. To obtain this tensor, sub¬ 
stitute from (60.4) into the equation c iikt = c mk and we 
obtain 

0* — v)(6 tk 8 n — d u 8 ik ) = 0. 

In this equation put i = k = l,j — l — 2 and we see that 
fi = v. Thus the isotropic tensor 

(61.1) c iik i = M" 8 kl + fx{d ik d)i + § (l 8j k ) 
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satisfies the required symmetry. The generalised Hooke’s 
law (59.1) for a homogeneous, isotropic body now 
becomes 

E (i = [AM« + bi&ikdii + £«£«)]«*! 
where A and ft are constants. This equation simplifies to 

(61.2) Eff = Xddjf -f 2/z£ w > 

where 9 is the dilatation defined by (56.5). Contracting 
i and j, we obtain the equation 


(61.8) 0 = E i( = (8A + 2/*)0 


which connects the two invariants 6 and 0. 

We now solve equations (61.2) for the strain tensor 
in terms of the stress tensor and in virtue of (61.8) we 
obtain 


(6l4) e “ = -viSITV) a «+^ £ «- 


It is usual to work in terms of Young’s modulus E and 
Poisson’s ratio a defined by 


(61.5) 


P _ /*(3A + 2/Q __ A 

a +/» ’ a ~2wnry 


Equation (61.4) then becomes 


(61.6) e it — -g { o0d ti -f* (1 4" o)E it }. 


We can obtain stress compatibility equations by sub¬ 
stituting this expression for e {j in (56.6). 

Sometimes we require the equations of equilibrium in 
terms of the displacement vector u t . To obtain these we 
substitute from (61.2) in (58.1) and obtain 

Ff + Xdidij + 2 fie iti = 0 
which becomes by (56.2) 

F t -f Xd ti + fi(u f ft + t if w ) = 0. 
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But 0 = e i{ = J(«<| + u i{ ) = u i{ . Consequently 
0 j = u i {i , and the atiove equations take the form 

(61.7) F, + (A + n)Q ti + = 0 

where P* is the Laplacian operator. Now that we have 
reverted to displacements u { we do not require any com¬ 
patibility equations. 

Ex. Show that if an isotropic, homogeneous body is in equilibrium 
under no forces, then P0 = 0. 

§ 62 . Curvilinear coordinates 
Many problems in elasticity can be investigated more 
conveniently by means of curvilinear coordinates x i in 
which the line-element da is given by da 2 = g {j dx i dx i . In 
this case we must return to the distinction between 
contravariance and covariance. 

The coordinate system serves only to describe the 
actual strains and stresses. The laws of elasticity are 
themselves independent of the coordinates. Therefore 
these laws can be formulated as tensor equations. We 
recall (section 9) that if a tensor is zero in any one co¬ 
ordinate system, it is zero in all coordinate systems. 
Consequently if we write down tensor equations, which 
reduce in cartesian coordinates to the results already 
established, then these tensor equations express the 
theory with respect to any curvilinear system. Ac¬ 
cordingly we immediately verify that the laws of elasticity 
are represented by the following tensor equations:- 

e = e if l l V, e u = £(«<.# + « M ), 6 = g«e w , 

e U,kl + e kl ,U e ik,il e il,ik = ®» 

T { = Eji n*, E if = Ejj, F* -f* E 1 ^ — 0, 

E{j = c {ikl eF l , c i} j.j = c ftk i = 

* It can also be proved that c kW = c iJtt . See I. S. Sokolnikoff, 
Mathematical Theory of Elasticity, section 26. 
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If the body is homogeneous and isotropic, we have in 
addition 

c Uhl = tyogkl + MSkiSli 4- £kigli)> 

E {i = Wg (i + 2 & = g u E ti = (8A + 2 fi)6, 

e u — {— oQ&a + (I 4* o)E it } 

F t + (* + n)B, t 4* = 0. 

In all these equations l* is a unit vector specifying the 
direction of the extension e , n* is the unit vector normal 
to the small element AS, and commas once more denote 
covariant differentiation. 

It is important to note that the components of the 
vectors u { , T { and F t may not possess a physical signi¬ 
ficance of the correct dimensionality. For example, the 
exercise of section 5 shows us that the second component 
of the acceleration vector in polar coordinates is an 
angular acceleration. It will suffice to discuss the force 
vector F it whose components in a cartesian system 

are F { = F l (say). Then F t = F. and F* = ~ F* 

are the covariant and contravariant components in the 
x* system, where we have put x { = y t the cartesian 
variables. The component of F*, which is the physical 
force vector, at any point in the direction of the 
unit vector l* = l ( is the invariant F*l { = FHf = 
g ii F i l i . The contravariant unit vector in the direction 
of the x 1 axis is dJ/Vin* Therefore the actual physical 
compone nts o f force along the x 1 coordinate curve are 
SuF'&ilV&I^gaF'IVri T- In the exercise of section 5, 
we have g u = 1, g ia = 0, g aa = r a and so the physical 
components of acceleration in polar coordinates are 
d* f /ddy , d*0 , n dr dd _. _ 

dfi ~ j and T ~di* + 2 dtlt’ whlch re P resent the 

radial and transverse accelerations respectively. 
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Similarly the tensors e u and E {j have no direct physical 
significance. Consider the strain tensor e i} which is 
connected with the cartesian strain tensor e kl by 

Qgk 

e ti — -~—t ; e k i. The physical component of the strain 

ox* ox* 

tensor associated with the directions of the unit vectors 
l* and m { at a point may be defined to be the invariant 
e t iVmf. 

As an example, let us discuss cylindrical coordinates 
r, 0 and z given by 

y 1 = r cos 0, y 2 = r sin 0, y s = z. 

The non-zero components of the fundamental tensor are 

Six ~ f>22 = r8 > gs8 = 1 

from which we deduce that 

S u — 1» g 22 =^, g 33 = 1, g« = 0ifi^/. 

So the only surviving Christoffel symbols of the second 
kind are 



The actual physical components (u a , up, Uy) of the dis¬ 
placement vector u t are 

1 

«« = «!> Up = — tt a , Uy = U a . 


Further calculation will show that the components of 
the strain tensor are 




CARTESIAN TENSORS - ELASTICITY 


103 


§ 63 


whilst the physical components associated with the direc¬ 
tions of the coordinate curves are 



The dilatation is 


0 = = 


du a du£ dtiy 
dr r 90 ~Bz 



Ex. Find the physical components associated with the directions 
of the coordinate curves of the strain tensor in terms of the 
physical components of the displacement vector when the 
coordinates are spherical polars. 


§ 68. Mechanics of continuous matter 
Let us discuss the motion of a continuous medium by 
the Eulerian method. Instead of following the path traced 
out by a particular particle, we focus our attention on a 
definite point P of the medium whose coordinates referred 
to a cartesian coordinate system are y { . Let us denote 
by tif the velocity vector of that particle which happens 
to be at P at time t. Then u, is a function of y t and t. 
After a further interval of time At has elapsed, the particle 
which was originally at P is now at the point y { + u ( At 
with velocity u { + Au ( . Hence u { + Au t is the function 
Uf at the point y { -f- u ( At and at time t -f- At. On expand¬ 
ing by Taylor’s theorem we have 

Mi -|- Ani -j- UjAtu { * -f* At • 

• dt 

The acceleration vector f { at P is the limit of AuJAt as 
At tends to zero. We therefore deduce that the acceleration 
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vector is given by 

( 68 . 1 ) f { = Ui u itt +^. 

Now consider a volume V of the medium, bounded by 
the surfaced.The massAf contained inside V is M=$ v pdV, 
where the density p is a function of y { and t. Thus the 

dl& C d 

rate of increase of mass is = I dV. Let n { denote 

at Jy Ot 

the direction cosines of the external normal to the small 
element of surface AS. Then the rate of mass flow out¬ 
wards across that element AS is pn { u { AS. Hence the rate 
of increase of mass is also expressed by the surface 
integral —J ' s pn ( u { dS. Consequently we have 

J pn { u { dS + J ^ dV = 0. 

We apply Gauss’s theorem to the surface integral and the 
result is 

J K (^).,dr+J r |« iv-o. 

This equation is an identity and so we derive the equa¬ 
tion of continuity 

( 68 . 2 ) ( pu t ),* + Tt =0 ‘ 

Physically this equation expresses the principle of con¬ 
servation of mass. 

In section 58 we discussed the equilibrium of a con¬ 
tinuous medium. The equations of motion can be deter¬ 
mined in the same way. If we replace F } by F f — pf } 
these equations become 

(68.8) pjj — Ff -j- E ii { . 

We now substitute for f t from (68.1) and obtain 
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ditt 

pu iUi ' t + p -l = F, + E liti . 

This equation can be written, in virtue of (63.2), in the 
form 

(68.4) [pUfUj — £</] i4 + ^ (pUf) = F f . 

(68.2) and (68.4) constitute the equations of motion of a 
continuous medium. 

When the coordinate system is curvilinear the equations 
(68.2) and (68.4) may be expressed in the tensor forms 

(68.5) (pu*)'i + \ = 0, 

(68.6) ( pu*u* — E”) it + A {put) = FK 


Solutions 


p. 88. Ex. The only non-vanishing components of a a are 
a ia “ — I* a ai = + 1» Oaa — 

^Ua 1 


p. 108. Ex. e M 


v ; 

1 lu. 


_ <>«£ «a . 
T 10 + T ’ 

u a cot 0 

H-+- 

r r 


r sin 0 

,* y =±p_^ + ±^7_22i® V 

2 \r sin 0 ’dtp r W r 

eya = + 

2 \r sin 0 by? 'br r J 

P 2 \r 7)0 T J>r r j 
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THEORY OF RELATIVITY 

§ 64. Special theory 

In flngg i rad mechanics, the position of a point in space 
at which an event occurs can be determined by its three 
space coordinates x\ x\ a 3 referred to some rectangular 
cartesian system. Also an observer can measure the time t 
at which the event takes place by means of a clock. An 
event is then fixed in both space and time by the system 
S which is comprised of the four numbers x 1 , ® 2 , jb 3 and t. 

Einstein examined the concept 'simultaneity* and came 
to the conclusion that ‘simultaneous events at different 
points’ has no meaning without further qualification. 
Continuing his study of fundamental ideas, Einstein 
arrived at the Special Theory of Relativity which he 
based on the two principles: (1) it is impossible to detect 
the unaccelerated translatory motion of a system through 
space, (2) the velocity c of a ray of light is a constant which 
does not depend on the relative velocity of its source and 
the observer. 

Let us now consider two systems 8 and 8 which coin¬ 
cide at the time t = 0, such that 8 moves with constant 
velocity V along the x 1 axis of the S system. Then the 
Lorentz Transformation, which can be deduced from 
the two principles of special relativity, connects the 
space coordinates and the time of both systems by the 
equations 

x 1 = ft(x l — Vt), x* = x\ z 3 = a?, i = (i(t — V&fc 2 ), 
where /S = (1—J^/c 2 ) - #. We easily verify that 
— (d* 1 ) 2 — (d® 2 ) 2 — (d® 3 ) 2 + c 2 (d*) 2 
= — (d® 1 ) 2 — (dm 2 ) 2 — (d® 3 ) 2 + c 3 (d*) 2 . 
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The invariance of this equation with respect to Lorentz 
transformations suggests that the Minkowski space 
defined by the metric 

(64.1) da 8 = — (<&>»)“ — (d® 8 ) 8 — (d® 8 ) 8 + c 8 (dc*) 8 

where we have written ®* = t, is appropriate for the 
geometrical discussion of special relativity. We denote the 
line-element of this four-dimensional space by da (not 
by ds) in order to emphasise that da is not the physical 
distance between two neighbouring points. 

The Minkowski space is flat and its signature, which 
equals the excess of the number of positive terms over 
the number of negative terms in its metric, is — 2. It 
is well-known* that there is no real transformation of 
coordinates which will reduce (64.1) to the metric of a 
four-dimensional Euclidean space, whose signature is 4. 
Thus the geometry of Minkowski space differs in many, 
respects from Euclidean geometry; for example, there 
exist real null curves (see (15.2)) and real null-geodesics. 

In this chapter Latin indices will have the range 1 to 8 
whilst Greek indices will range from 1 to 4. The velocity 
« of a particle, which is at the point x* has the components 
u* = da { /dt referred to the system S. It follows from 
(64.1) that 



The four-dimensional Minkowski momentum vector 

is defined by m 0 c , where m 0 is a constant. The special 

** da 1 

theory identifies the fourth component with the 

mass m of the moving particle. In virtue of (64.2) we 
have 

(64.8) m = m 0 (l — *. 

* W. L. Ferrar, Algebra, p. 154. 
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The constant m 0 is the mass when tt = 0 and so it is 
called the rest-mass of the particle. The mass m, which 
clearly increases with the velocity, is called the rela¬ 
tivistic mass of the particle. The components 


m,c 


dx*_ 

da 


= m 0 c 


das* dt _ das* 
1ido~ m !t 


and are evident generalizations of the Newtonian mo¬ 
mentum vector. 

We define the four-dimensional Minkowski force 
vector F a by 


(64.4) 


9 d 2 ®“ » dl d®“\ /, tt 2 \ - Hd/ dx°\ 

Fa - m ^~do s ’- c 5;\ mo do) - ( 1— dt[ m nr) 

d / d&*\ 

The Newtonian force vector is X* = — (m —I and sc 

We obtain by expansion from (64.8) that 


WIC 2 = 77*0C a + — OTott 8 + ... 


and so the special theory identifies the energy E as¬ 
sociated with a particle by the equation E — me 2 . 
Therefore 

L «V H dm 1 / « 2 \-i4 dE 

The motion of a particle which moves under the action 
of some force system can be represented in Minkowski 
space by a curve, called the world-line of the particle. 
If no forces act on the particle, we see from (64.4) that 
d 2 ®“/do 2 = 0. Thus the world-line of a free particle is a 
geodesic of the Minkowski space. 

The velocity of a light ray is the constant c, and so 
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we see from (64.1) that for such a ray da = 0. Accor¬ 
dingly the world-line of a light ray is a null-geodesic of 
the Minkowski space. 

In order to discuss the mechanics of a continuous 
medium, we introduce the symmetrical four-dimensional 
energy-momentum tensor T*0 defined by 

= T H = pu*U> — E*U T li = T“ = pu<; T 44 = p, 

where p is the density and E** is the Cartesian stress 
tensor defined in section 67. Then the special theory 
generalises (68.5) and (68.6), which are the equations of 
motion of a continuous medium into 

(64.5) T* = Ffi. 

If we change to spherical polar coordinates r, 6 and ip, 
the metric of Minkowski space becomes 

(64.6) da* — — dr* — r*dd 2 — r 2 sin* Qdip* c 2 di*. 


§ 65. Maxwell’s Equations 
The classical theory of electrodynamics*, according 
to Lorentz, is specified by the electric potential q> which 
is a scalar and the magnetic potential A { which is a 
vector. The electric field strength vector E { and the 
magnetic field strength vector H { are derived from 
these potentials by the equations 

j? , 1 dA { 

E { -— grad <p -=-*, 

c ot 


H { = curl 


* M. Abra ham —R. Becker, Electricity and Magnetism* 

C. A* Coulson, Waves, chap. VII. 

D. E. Rutherford, Vector Methods, p. 126. 
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Using electrostatic units, Maxwell’s equations are 


(65.1) 


div E t — 4otp, 
div H { — 0, 


0011 b ‘ + 71 t‘“ 0 ' 

. _ 1 dE { 4n . 


where j t is the current density vector and p is the 
charge density. 

In Minkowski space, with the metric (64.1), let us form 
the four-dimensional potential vector 0 a and the four- 
dimensional current density vector J° defined respec¬ 
tively by 

0a = (— -di» — ^2» — 

— 0"i > 7 a > 7a» p)> 


with respect to a particular coordinate system. Next we 
introduce the skew-symmetrical tensor t] a p defined by 


Va0 — 0a J 0fi,a 


d0a 

dafi 


d0fi 

dx* 


and we immediately calculate that its non-va n is hin g 
components in the given coordinate system are 

Vts = Vaa = E i> Vai~ Via = E a> *?i2 = Vn = E a» 
Vu = ~Va = cE v Vai = — Vi2 = cE a> Vzi = ~Vi& = cE a- 

The non-vanishing contravariant components rfP may 
now be obtained and are 


V 

V 


23 =— r) sz = H x ; rf l — — ij 1 * = H 3 ; —jj 81 = ff 3 ; 

14 = _ ?? U = _^l. ^34 = _^43 = _^3. 


We now write Maxwell’s equations (65.1) in terms of 
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rj and J and the results are readily verified to be res¬ 
pectively 

dr]* 1 dr]** dr)** _ 4m 4 
dx 1 dx* da? ~ c * 

d Vsa , si , fyu n 
dx 1 "*■ da? + da? ~ ’ 

dy<t i dVii , dVv _ A 
da? ^ dx* * dx> ~ ’ 
drf 1 dt) i2 dt] i3 dr])* _ 4m J{ 

da? da? da? da? c 

The first and last of these equations combine together 
into the form 

( 65 . 2 ) 

whilst the remaining two are accounted for by the 
equations of the set 

( 65 . 8 ) rj affi7 + t) fiVia + r] yatfi = 0 

which do not vanish identically. 

We have accordingly written Maxwell’s equations in 
tensor form in Minkowski space. Thus they are invariant 
under the Lorentz group of transformations. 

§ 66. General theory 

We now turn to the General Theory of Relativity which 
was developed by Einstein in order to discuss gravitation. 
He postulated the principle of covariance, which 
asserts that the laws of physics must be independent of 
the space-time coordinates. This swept away the privi¬ 
leged role of the Lorentz transformation. As a result 
Minkowski space was replaced by the Riemannian V t with 
the general metric 

( 66 . 1 ) do* = gafidopda?. 

Einstein also introduced the principle of equivalence. 



112 


TENSOR CALCULUS 


§ 66 

which in essence states that the fundamental tensor 
g a fi can be chosen to account for the presence of a gravi¬ 
tational field. That is, g a $ depends on the distribution of 
matter and energy in physical space. 

Matter and energy can be specified by the energy- 
momentum tensor T a & which in the special theory satisfies 
the equation T^j ,f = F#. The only forces, namely those 
due to gravitation, are however already taken into 
account by the choice of the fundamental tensor g a p. 
We therefore ignore F$ and, in accordance with the prin¬ 
ciple of covariance, the energy-momentum tensor must 
now satisfy the equation T*jj = 0. We shall write this 
equation in the equivalent form a = 0 where 
T% — ifiyT a v is the mixed energy-momentum tensor. 
The problem now is to determine as a function of 
the g a ft and their derivatives up to the second order, 
bearing in mind that a = 0. We recall from (84.8) 
that Einstein’s tensor defined by 

(66.2) = g«r R ffr — 

satisfies the equation G*p a — 0. The equations of motion 
require a = 0, but very remarkably a = 0 is 
an identity in Riemannian geometry. This led Einstein to 
propose the relation 

(66.8) kT% + G* fi = 0. 

In effect these equations form the link between the 
physical energy-momentum tensor and the geome¬ 
trical tensor of the F 4 of general relativity. In order 
that Newton’s theory of gravitation can be deduced as 
a first approximation from Einstein’s theory, it was 
found necessary to choose k = inkjet where k is the 
usual gravitational constant 6.664 x 10“ 8 cm 3 , gm. -1 
sec. -8 . The value of c is 2.99796 x 10 10 cm. sec. -1 , and 
so k is 2.078 x 10 -48 cm. -1 gm. -1 sec. 8 in e.g.s. units. 
In the special theory, the world-lines of free particles 
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and of light rays are respectively the geodesics and the 
null-geodesics of Minkowski space. The principle of equi¬ 
valence demands that all particles be regarded as free 
particles when gravitation is the only force under con¬ 
sideration. Then it follows from the principle of co- 
variance that the world-line of a particle under the action 
of gravitational forces is a geodesic of the F 4 with the 
metric (66.1). Similarly the world-line of a light ray is a 
null-geodesic. 


§ 67. Spherically symmetrical metric 
General relativity discusses several important problems 
in which the coordinate system r, 0, ip and t is such that 
the metric takes the form 

(67.1) da 2 = — e A dr 2 — r a dO 2 — r a sin 2 6 dtp 2 + c 2 e v dt 2 , 

where A and v are functions of r. A metric of this type 
is said to be spherically symmetrical. It is a generalisa¬ 
tion of the special relativity metric (64.6), which is ex¬ 
pressed in spherical polars. The coefficients of dr 2 and dt a 
have been selected as exponentials in order to ensure 
that the signature of da 2 is — 2. Let us write x 1 = r, 
x 2 = 6, x 2 = y> and x* = ct. Then the non-zero components 
of the fundamental tensor are 


gu = — e\ g 22 = — r 2 , g 33 = — r 2 sin 2 0, g M = e v . 
The determinant g becomes 

g = — r* sin 8 0 e x+p 


and hence the non-zero components of the conjugate 
symmetric fundamental tensor are 


»ll = _ e -A 



r>33 


1 

r 2 sin 2 0 ’ 


g« = e~ v . 


A brief calculation shows that the only non-vanishing 
Christoffel symbols of the second kind are 
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where a dash denotes differentiation with respect to r. We 
now evaluate the components of the Ricci tensor by 
means of (88.2) and the results are 


Rn - - j A'- *AV + \v" + 

(67.8) ^ 2a = cose<j8 0 ®ss = — 1 + $rv'}, 

R u = — \v"—jv’ — Jv'a}, 

R a fi — 0 for a # 0. 

A further calculation gives us the curvature invariant 

(67.4) . 

On substituting this in (66.2) we obtain the components 
of the Einstein tensor for the spherically symmetrical 
metric (67.1) 



Ex. Find the necessary and sufficient conditions that a space with 
a spherically symmetrical metric be an Einstein space. 
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§ 68. Schwarzschild metric 

We now seek the spherically symmetrical metric 

(67.1) consistent with the existence of one gravitating 
point particle situated at the origin, and surrounded by 
empty space. When the origin itself is excluded from our 
discussion, the energy-momentum tensor is zero at 
all points. It follows from (66.8) that G*p = 0, which 
yields in virtue of (67.6) the equations 

( 68 . 1 ) + + 

(68.2) «r»{- ll ’-iiV + | V" +1,' + iy) = 6, 

(68.8) -I = 

The solution of (68.8) is readily found to be e-*=1— 2m/(c 2 r), 
where the constant of integration m introduced in this 
way can be identified physically with the rest mass of the 
gravitating particle. Further, we subtract (68.8) from 
(68.1) and obtain e~*(X' + v')(r = 0. That is X' v' = 0, 
from which we have X 4- v = k, where As is a constant. 
Thus e v = e k {l — 2m/(c 2 r)}, We can now verify that 
equation (68.2) is satisfied identically. However, at large 
distances from the gravitating point, the metric should 
approximate to the metric (64.6) of special relativity. 
Therefore, we must select k = 0. We thus obtain the 
Schwarzschild metric 


(68.4) 


do* = 



dr 8 — r 8 dQ 2 — r a sin 8 Odyfl 



Ex. Show that a space with Schwarzschild’s metric is an Einstein 
space, but not a space of constant curvature. 
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§ 69. Planetary motion 

Let us investigate the motion of a planet in the gravita¬ 
tional field of the sun. The sun will be selected as a 
gravitating particle and the planet as a free particle 
whose mass is so small that it does not affect the metric, 
and whose world-line is then a geodesic in the V t with the 
Schwarzschild metric (68.4). The geodesics are deter¬ 
mined by the four equations (26.4) in which we now, 
of course, replace 8 by a. We shall omit one of these 
equations, in practice the most formidable one involving 
tPr/do 8 , and replace it by (27.1) which is of the first order 
and is satisfied along the geodesics. When this is done, 
we no longer require the Christoffel symbols of the type 

. The remaining non-vanishing symbols of the second 

kind can be calculated from (67.2) and may be found 
to be 



12 1 _ jl_ f4 1 _ m / 2m\-i 

Iwj “7* \l8j7’ \l4j \ 7r/ ' 

t 2 \ ( Q \ 

881 = — sin 0 cos l 2 8| = COt 


Hence the four equations of the geodesics are 
... _. d 9 0 2 dr dd . (dyV 

(691) ^ + 7*S- s “‘ 9om6 (s) =*• 


i 2T + l^ d JP + 2mt e d ^f = o. 

da 2 r da da da da 
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We may assume that the planet moves initially in the 
plane 0 = Jt/2. That is, dd/do and cos 0 are both initially 
zero. Then (69.1)tells us that d a 0/d<7 a is also zero. Repeated 
differentiation of this equation shows that d l d/do ( vanishes 
at t = 0 for all i. Hence 6 = jt/2 permanently, and the 
above equations simplify to 


(69.2) 

d 2 tp 2 dr dip 
da 2 r da da 


(69.8) 

d?t 2m / 2m\-idr di _ 

d<y a c a r \ c a r/ da da 


(69.4) - 


2m\ / di \ 8 
c a r/ \da) 


We can immediately integrate (69.2) and (69.8) and the 
results are 


(69.5) 



where h and k are constants. On eliminating t and or 
from (69.4) and (69.5) we obtain 



Now substitute r = 1/t* and differentiate the equation so 
obtained with respect to y> and the result is 


(69.6) 


d?u m 


Zmu* 

~cT 


For the planets of our solar system, the term m/c a A 2 is 
much larger than 8mu 2 /c 2 . But when we neglect this latter 
term, we obtain Newton’s equation for the motion of a 
planet. Thus the first approximation to the solution of 


(69.6) is the Newtonian solution u = {l+£cos(y—£)}, 


where e is the eccentricity of the elliptic orbit and £ is 
the longitude of perihelion. A second approximation 
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to the solution can then be obtained in the form 
u = J3/ja I 1 + e 008 (V—f— d£)}> where A£ = BmPf/c^h 2 . 

This means that the major axis of the elliptic orbit is 
slowly rotating about its focus (the sun). The increase 
in AS corresponding to a complete revolution y> = 2jc is 
thus Qmbifcth 2 . For the planet Mercury the advance of 
perihelion is calculated from this to be 42.9 seconds of 
arc per century. This agrees well with the observational 
figures of 48.5 seconds of arc per century. 


§ 70. Einstein’s universe 

Einstein was led by cosmological considerations to con* 
sider the universe with the metric 

(70.1) da 2 = —(1—r 2 /^ 2 ) -1 dr 2 — r 2 d0 2 —r a sin 2 6dip 2 -\-c 2 di 2 , 

where 01 is a constant. This metric is spherically symme¬ 
trical with e~*- = (1 — r 2 ffi 2 ) and v = 0. The Christoff el 
symbols of the second kind are readily obtained from 

(67.2) . 

Let us investigate the path of a ray of light in Ein¬ 
stein’s universe. The path must be a null-geodesic and 
so its equations are given by three of the four equations 

(26.4) taken together with = 0. That is, we 

have du du 


(70.2) 


d?6 2 dd dr 

•js H-j~j -sin 0 cos 0 

au a r au flu 


$)’=»• 


2 dip dr dip dd 

j a H-j 2 - j- + 2 cot 0 -p- — = 0, 

flu 2 r flu du du du 


dH 
du 2 


= 0 , 


—(1—r 2 /^ 2 )-i(^Y—r 2 ®—r 2 sin 2 0(^Y+c 2 (^y= 0, 

\duf \duj \du/ \duj 
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$ 70 


where u is some parameter. Again, following the argument 
of section 69, the equation (70.2) tells us that we can 
take 6 to have the permanent value ?r/2. With this choice, 
the remaining equations reduce to 


(70.8) 


(Pip 2 dip dr _ d a t 
dtp r du du ’ dtp 


(«••*) - —■*/***>“ (as)’ - r ' {ZJ+ c ‘ (s)‘“ °- 

We integrate (70,8) and obtain 


(70,5) 




where h and k are constants, and then eliminate t and u 
from these equations and (70.4). The result is 



The solution of this equation is 


1 1 jfc a 

(70.6) - = — cos* {ip —£) + -^ -sin*(v> —I) 

where f is a constant. We immediately see that r regains 
its initial value when ip is increased by n and that r is 
never infinite for any value of ip. Thus all the null* 
geodesics of Einstein’s universe, that is the light rays, 

dt k 

are closed curves. We see from (70.5) that r a . 

Hence the time taken for a light ray to make a complete 
circuit is given by 

k r *» { 1 c s k 2 1 _1 

T=jj o {^5 cosa (V — f) + -p-sin*(v> — f)J dip. 
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Owing to the periodicity of ip we have 

-ini— 

4ft f"/8 t 1 , . c 2 ft a . 

= tJ 0 [m> cos *v 


sin 2 ip\ dip 


and on carrying out the integration we find thatT=2jr^/c. 

Ex. 1. Show that the Einstein universe is neither an Einstein 
space nor a space of constant curvature. 

Ex. 2. Prove that the curvature invariant of Einstein’s universe 
is R = 6 \3i\ 


§ 71. De Sitter’s universe 
Other cosmological considerations suggested to De 
Sitter that the universe could be described by the metric 

(71.1) do 2 = — (1 — r 2 /^ 2 )" 1 dr 2 —rW—r a sin 2 0 dip 2 

+ c 2 (l—r 2 /^ a )df 2 . 

This metric has also spherical symmetry but the con¬ 
stant 0t has not the same value as the corresponding 
constant of the Einstein universe. 

The paths of light rays are the null-geodesics given by 
the equations 


d a 0 ,2 dr dd . n n 

-j-j H--j- j -sm 0 cos 0 

dti 2 r du du 


(M* _ o 

\du) 


+ + o. 

< 71 - 8 > £-£o— 


(71.4) — (1 


a /^ 2 ) -1 (tV—ir 2 r 2 sin 2 0 feV 

\au/ \au/ \aw/ 

+ C 2(l_ r 2/^2)^*j 2 = 0. 
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5 71 

Again we can choose 0 = jtj2 permanently, and on in¬ 
tegration the equations (71.2) and (71.8) become 

r ‘Tu = h - | = 

We now eliminate t and u from these equations and 
(71.4). The result is 



where a 8 = c 8 ft 8 /A 8 + l/^ 8 . This equation can be im¬ 
mediately integrated and yields 

(71.5) — = a cos (ip — £), 

T 

where £ is a constant. These trajectories correspond to 
straight lines and are not closed, since r becomes infinite 
when y> — £ = n/2. 

Ex. Show that the De Sitter universe is an Einstein space with 
constant curvature 12/&*. 

Solution 

p. 114. Ex. e - * = 1 + ar * + b/r ; v = k — A, where a, b and k 
are constants. 
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